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PREFACE 

This workbook presents some computational techniques currently used by scien- 
tists working with atmospheric dispersion problems. Because the basic working equa- 
tions are general, their application to specific problems usually requires special care 
and judgment; such considerations are illustrated by 26 example problems. This 
workbook is intended as an aid to meteorologists and air pollution scientists who are 
required to estimate atmospheric concentrations of contaminants from various types 
of sources. It is not intended as a complete do-it-yourself manual for atmospheric 
dispersion estimates; all of the numerous complications that arise in making best esti- 
mates of dispersion cannot be so easily resolved. Awareness of the possible complex- 
ities can enable the user to appreciate the validity of his "first approximations" and 
to realize when the services of a professional air pollution meteorologist are required. 
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ABSTRACT 

This workbook presents methods of practical application of the binormal con- 
tinuous plume dispersion model to estimate concentrations of air pollutants. Estimates 
of dispersion are those of Pasquill as restated by Gifford. Emphasis is on the estima- 
tion of concentrations from continuous sources for sampling times up:, to 1 hour.. Some 
of the topics discussed are determination of effective height of emission, extension of 
concentration estimates to longer sampling intervals, inversion break-up fumigation 
concentrations, and concentrations from area, line, and multiple sources. Twenty-six 
example problems and their solutions are given. Some graphical aids to computation 
are included. 
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Chapter 1 INTRODUCTION 



During recent years methods of estimating at- 
mospheric dispersion have undergone considerable 
revision, primarily due to results of experimental 
measurements. In most dispersion problems the 
relevant atmospheric layer is that nearest the 
ground, varying in thickness from several hundred 
to a few thousand meters. Variations in both 
thermal and mechanical turbulence and in wind 
velocity are greatest in the layer in contact with 
the surface. Turbulence induced by buoyancy forces 
in the atmosphere is closely related to the vertical 



temperature structure. When temperature decreases 
with height at a rate higher than 5.4F per 1000 ft 
(1C per 100 meters), the atmosphere is in un- 
stable equilibrium and vertical motions are en- 
hanced. When temperature decreases at a lower 
rate or increases with height (inversion), vertical 
motions are damped or reduced. Examples of typ- 
ical variations in temperature and wind speed with 
height for daytime and nighttime conditions are 
illustrated in Figure 1-1. 
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Figure 1-1. Examples of variation of temperature and wind speed with height (after Smith, 1963). 



The transfer of momentum upward or down- 
ward in the atmosphere is also related to stability; 
when the atmosphere is unstable, usually in the 
daytime, upward motions transfer the momentum 
"deficiency" due to eddy friction losses near the 
earth's surface through a relatively deep layer, 
causing the wind speed to increase more slowly 
with height than at night (except in the lowest few 
meters). In addition to thermal turbulence, rough- 
ness elements on the ground engender mechanical 
turbulence, which affects both the dispersion of 
material in the atmosphere and the wind profile 
(variation of wind with height). Examples of these 
effects on the resulting wind profile are shown in 
Figure 1-2. 



As wind speed increases, the effluent from a 
continuous source is introduced into a greater vol- 
ume of air per unit time interval. In addition to 
this dilution by wind speed, the spreading of the 
material (normal to the mean direction of trans- 
port) by turbulence is a major factor in the dis- 
persion process. 

The procedures presented here to estimate at- 
mospheric dispersion are applicable when mean wind 
speed and direction can be determined, but meas- 
urements of turbulence, such as the standard de- 
viation of wind direction fluctuations, are not avail- 
able. If such measurements are at hand, techniques 
such as those outlined by Pasquill (1961) are likely 
to give more accurate results. The diffusion param- 



eters presented here are most applicable to ground- 
level or low-level releases (from the surface to about 
20 meters), although they are commonly applied at 
higher elevations without full experimental valida- 
tion. It is assumed that stability is the same 
throughout the diffusing layer, and no turbulent 
transfer occurs through layers of dissimilar stability 
characteristics. Because mean values for wind direc- 
tions and speeds are required, neither the variation 
of wind speed nor the variation of wind direction 
with height in the mixing layer are taken into ac- 
count. This usually is not a problem in neutral or 
unstable (e.g., daytime) situations, but can cause 
over-estimations of downwind concentrations in 
stable conditions. 
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ATMOSPHERIC DISPERSION ESTIMATES 



Chapter 2 BACKGROUND 



For a number of years estimates of concentra- 
tions were calculated either from the equations of 
Sutton (1932) with the atmospheric dispersion 
parameters C y , C, and n, or from the equations of 
Bosanquet (1936) with the dispersion parameters 
p and q. 

Hay and Pasquill (1957) have presented experi- 
mental evidence that the vertical distribution of 
spreading particles from an elevated point is re- 
lated to the standard deviation of the wind eleva- 
tion angle, <r,. ; , at the point of release. Cramer ( 1957) 
derived a diffusion equation incorporating standard 
deviations of Gaussian distributions: _,- for the 
distribution of material in the plume across wind 
in the horizontal, and ^ for the vertical distribution 
of material in the plume. (See Appendix 2 for prop- 
erties of Gaussian distributions.) These statistics 
were related to the standard deviations of azimuth 
angle, <r A , and elevation angle, <r,. ; , calculated from 
wind measurements made with a bi-directional 
wind vane (bivane). Values for diffusion param- 
eters based on field diffusion tests were suggested 
by Cramer, et al. (1958) (and also in Cramer 1959a 
and 1959b). Hay and Pasquill (1959) also pre- 
sented a method for deriving the spread of pollut- 
ants from records of wind fluctuation. Pasquill 
(1961) has further proposed a method for esti- 
mating diffusion when such detailed wind data are 
not available. This method expresses the height 
and angular spread of a diffusing plume in terms of 
more commonly observed weather parameters. Sug- 
gested curves of height and angular spread as a 
function of distance downwind were given for sev- 
eral "stability" classes. Gifford (1961) converted 
Pasquill's values of angular spread and height into 
standard deviations of plume concentration distri- 
bution, ff y and (T K . Pasquill's method, with Giflord's 
conversion incorporated, is used in this workbook 
(see Chapter 3) for diffusion estimates. 

Advantages of this system are that (1) only two 
dispersion parameters are required and (2) results 
of most diffusion experiments are now being re- 
ported in terms of the standard deviations of plume 
spread. More field dispersion experiments are being 
conducted and will be conducted under conditions 
of varying surface roughness and atmospheric sta- 
bility. If the dispersion parameters from a specific 
experiment are considered to be more representative 



than those suggested in this workbook, the param- 
eter values can be used with the equations given 
here. 
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Chapter 3 ESTIMATES OF ATMOSPHERIC DISPERSION 



is chapter outlines the basic procedures to 
sd in making dispersion estimates as sug- 
by Paaquill (1961) and modified by Gifford 

DINATE SYSTEM 

the system considered here the origin is at 
I level at or beneath the point of emission, 
ie x-axis extending horizontally in the direc- 
1 the mean wind. The y-axis is in the hori- 
plane perpendicular to the x-axis, and the 
extends vertically. The plume travels along 
illel to the x-axis. Figure 3-1 illustrates the 
late system. 

SIGN EQUATIONS 

3 concentration, x, of gas or aerosols (parti- 
is than about 20 microns diameter) at x,y,z 
continuous source with an effective emission 
H, is given by equation 3.1. The notation 
o depict this concentration is x (x,y,z;H). 
he height of the plume centerline when it 



becomes essentially level, and is the sum of the 
physical stack height, h, and the plume rise, AH. 
The following assumptions are made: the plume 
spread has a Gaussian distribution (see Appendix 
2) in both the horizontal and vertical planes, with 
standard deviations of plume concentration distri- 
bution in the horizontal and vertical of o- y ^and a zs 
respectively; the mean wind speed affecting the 
plume is u; the uniform emission rate of pollutants 
is Q; and total reflection of the plume takes place 
at the earth's surface, i.e., there is no deposition 
or reaction at the surface (see problem 9). 



X (x,y,z;H) 



exp 



z+H 



Q 



exp' 



z-H 



a 

1+ exp 1 
1* 



1 
2 

(3.1) 



'Note: exp a/b = e- n />> where e Is the base of natural logarithms 
and is approximately equal to 2.7183. 
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Figure 3-1. Coordinate system showing Gaussian distributions in the horizontal and vertical. 
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Any consistent set of units may be used. The most 
common is; 

\ (g m" 1 ) or, for radioactivity (curies m~ 3 ) 

Q (gsec~ l ) or (curies sec' 1 ) 

u ( m sec~ ! ) 

IT,, a tt H,x,y, and z (m) 

This equation is the same as equation (8.35) p. 293 
of Sutton (1953) when u's are substituted for But- 
ton's parameters through equations like (8.27) p. 
286. For evaluations of the exponentials found in 
Eq. (3.1) and those that follow, see Appendix 3. 
x is a mean over the same time interval as the time 
interval for which the a's and u are representative. 
The values of both a y and tr t are evaluated in terms 
of the downwind distance, x. 

Eq. (3.1) is valid where diffusion in the direc- 
tion of the plume travel can be neglected, that is, 
no diffusion in the x direction. 

This may be assumed if the release is continuous 
or if the duration of release is equal to or greater 
than the travel time (x/u) from the source to the 
location of interest. 

For concentrations calculated at ground level, 
i.e., z = 0, (see problem 3) the equation simplifies 
to: 



X (x,y,0;H) = 



Q 



exp 



(3.2) 

Where the concentration is to be calculated 
along the centerline of the plume (y = 0), (see 
problem 2) further simplification results: 



X <x,0,0;H) 






LJJ'JLVI 

I 2 UJ ](3.3) 



_ For a ground-level source with no effective plume 
rise (H = 0), (see problem 1): 



X (x,0,0;0) 



Q 



U 



(3.4) 



EFFECTS OF STABILITY 



The values of <r y and CTZ vary with the turbulent 
structure of the atmosphere, height above the sur- 
face, surface roughness, sampling time over which 
the concentration is to be estimated, wind speed 
and distance from the source. For the parameter 
values given here, the sampling time is assumed to 
be about 10 minutes, the height to be the lowest 
several hundred meters of the atmosphere, and 
the surface to be relatively open country. The 
turbulent structure of the atmosphere and wind 
speed are considered in the stability classes pre- 



sented, and the effect of distance from the soums IH 
considered in the graphs determining the parameter 
values. Values for <r y and <r B are estimated from the 
stability of the atmosphere, which is in turn esti- 
mated from the wind speed at a height of about 
10 meters and, during the day, the incoming solar 
radiation or, during the night, the cloud cover (Pan- 
quill, 1961). Stability categories (in six classen) 
are given in Table 3-1. Class A is the most un- 
stable, class F the most stable class considered horo. 
Night refers to the period from 1 hour before aunHot 
to 1 hour after sunrise. Note that the neutral 
class, D, can be assumed for overcast condition** 
during day or night, regardless of wind speed. 

Table 3-1 KEY TO STABILITY CATEGORIES 
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The neutral class, D, should be assumed for overcast conditions durlnfi 
day or night. 

"Strong" incoming solar radiation corresponds 
to a solar altitude greater than 60 with clear skies; 
"slight" insolation corresponds to a solar altitude; 
from 15 to 35 with clear sides. Table 170, Solar 
Altitude and Azimuth, in the Smithsonian Mete- 
orological Tables (List, 1951) can be used in deter- 
mining the solar altitude. Cloudiness will decrease! 
incoming solar radiation and should be considered 
along with solar altitude in determining solar radia- 
tion. Incoming radiation that would be strong 
with clear skies can be expected to be reduced to 
moderate with broken (% to % cloud cover) mid- 
dle clouds and to slight with broken low cloudH 
An objective system of classifying stability from 
hourly meteorological observations based on the 
above method has been suggested (Turner, 1961). 

These methods will give representative indica- 
tions of stability over open country or rural areas 
but are less reliable for urban areas. This differ- 
ence is due primarily to the influence of the city's 
larger surface roughness and heat island effects 

Zl ^ S ff tabiHty * over urban areas .Tte 
greatest difference occurs on calm dear nights- on 

s S ta C ble m h g \ tS C ndit r 8 Ver mral areas are ^y 
stab e, but over urban areas they are slightly un- 
stable or near neutral to a height several times the 
average building height, with a stable lajwTboro 
(Duckworth and Sandberg, 1954; DeMarrais'igeT) 
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Some preliminary results of a dispersion experi- 
ment in St. Louis (Pooler, 1965) showed that the 
dispersion over the city during the daytime behaved 
somewhat like types B and C; for one night, experi- 
ment o-j. varied with distance between types D and E. 

ESTIMATION OF VERTICAL AND 
HORIZONTAL DISPERSION 

Having determined the stability class from 
Table 3-1, one can evaluate the estimates of o-,. and 
<r K as a function of downwind distance from the 
source, x, using Figures 3-2 and 3-3. These values 
of CTV and try. are representative for a sampling time 
of about 10 minutes. For estimation of concentra- 
tions for longer time periods see Chapter 5. Figures 
3-2 and 3-3 apply strictly only to open level country 
and probably underestimate the plume dispersion 
potential from low-level sources in built-up areas. 
Although the vertical spread may be less than the 
values for class F with very light winds on a clear 
night, quantitative estimates of concentrations are 
nearly impossible for this condition. With very light 
winds on a clear night for ground-level sources free 
of topographic influences, frequent shifts in wind 
direction usually occur which serve to spread the 
plume horizontally. For elevated sources under 
these extremely stable situations, significant con- 
centrations usually do not reach ground level until 
the stability changes. 

A stable layer existing above an unstable layer 
will have the effect of restricting the vertical diffu- 
sion. The dispersion computation can be modified 
for this situation by considering the height of the 
base of the stable layer, L. At a height 2.15 
above the plume centerline the concentration is one- 
tenth the plume centerline concentration at the same 
distance. When one-tenth the plume centerline 
concentration extends to the stable layer, at height 
L, it is reasonable to assume that the distribution 
starts being affected by the "lid." The following 
method is suggested to take care of this situation. 
Allow (r z to increase with distance to a value of 
L/2.15 or 0.47 L. At this distance x r ,, the plume is 
assumed to have a Gaussian distribution in the 
vertical. Assume that by the time the plume travels 
twice this far, 2 XL, the plume has become uniformly 
distributed between the earth's surface and the 
height L, i.e., concentration does not vary with 
height (see Figure 3-4). For the distances greater 
than 2 xi,, the concentration for any height between 
the ground and L can be calculated from: 



X (x,y,z;H) 



Q 



V2rr o-y Lu 



exp I 



(3.5) 



(see problem 6). Note that Eq. (3.5) assumes nor- 
mal or Gaussian distribution of the plume only in 
the horizontal plane. The same result can be ob- 
tained from the following equation where tr zlj is an 
effective dispersion parameter because V^TT L = 
2.5066 L and 0.8 TrL = 2.51 L. 



X (x,y,z;H) 



Q 



71- <r v 0-j.L U 



exp 



(3.6) 



for any z from to L 
for x >2 XT,; x r , is where o- K 



0.47 L 



for any z from to L 

for x > 2. v r.; x [: is where u- y . =- 0.47 L 
The value of <r K] , = 0.8 L 

EVALUATION OF WIND SPEED 

For the wind speed, u, a mean through the ver- 
tical extent of the plume should be used. This 
would be from the height H 2 v z through H -f- 
2 CT K . Of course, if 2 o- K is greater than H then the 
wind can be averaged from the ground to H -f- 2 a x . 
However, the "surface wind" value may be all that 
is available. The surface wind is most applicable 
to surface or low-level emissions, especially under 
stable conditions. 

PLOTS OF CONCENTRATIONS 
AGAINST DISTANCE 

To gain maximum insight into a diffusion prob- 
lem it is often desirable to plot centerline concen- 
trations against distance downwind. A convenient 
procedure is to determine the ground-level center- 
line concentrations for a number of downwind dis- 
tances and plot these values on log-log graph paper. 
By connecting the points, one may estimate con- 
centrations for intermediate downwind distances 
(see problem 6). 

ACCURACY OF ESTIMATES 

Because of a multitude of scientific and techni- 
cal limitations the diffusion computation method 
presented in this manual may provide best estimates 
but not infallible predictions. In the unstable and 
stable cases, severalfold errors in estimate of <T K 
can occur for the longer travel distances. In some 
cases the a* may be expected to be correct within a 
factor of 2, however. These are: (1) all stabilities 
for distance of travel out to a few hundred meters; 
(2) neutral to moderately unstable conditions for 
distances out to a few kilometers; and (3) unstable 
conditions in the lower 1000 meters of the atmos- 
phere with a marked inversion above for distances 
out to 10 km or more. Uncertainties in the esti- 
mates of tr y are in general less than those of <r z . 
The ground-level centerline concentrations for these 
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Figure 3-2. Horizontal dispersion coefficient as a function of downwind distance from the 



source. 
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Figure 3-3. Vertical dispersion coefficient as a function of downwind distance from the source. 
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Figure 3-4. Variations in concentration in the vertical beneath a more stable layer. 



three cases (where ^ can be expected to be within 
a factor of 2) should be correct within a factor of 3 
including errors in Vy and u. The relative confidence 
in the ,s (m decreasing order) is indicated by the 
heavy lines and dashed lines in Figures 3-2 and 3-3. 

Estimates of H, the effective height of the plume, 
may he in error because of uncertainties in the esti- 
mation of AH the plume rise. Also, for problems 
that require estimates of concentration at a specific 
point the difficulty of determining the mean" wind 
" time interval d consequently the 
X ' aXiS Can un- 



PLOTTING GROUND-LEVEL 
CONCENTRATION ISOPLETHS 

Often one wishes to determine the locatioiiH 
where concentrations equal or exceed a given mag- 
nitude. First, the axial position of the plume must 
be determined by the mean wind direction For 
plotting isopleths of ground-level concentrations 
the relationship between ground-level centorlino 
concentrations and ground-level off-axis concentra- 
tions can be used: 



.x (*,y,0;H) 
x (x,0,0;H) 



GRAPHS FOR ESTIMATES OF DIFFUSION J! J oor ? ina f e f a particular isopleth from the 

*' A-axis can be dprpvmmoi-i ^v- nn i. .i_. 
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times wind speed (v u/Q) aeainst down 
wind distances for various effective heifhts of * ml" 
sion and limits to the vertical mixing foi eachT4 
bihty class (1 figure for each stability) Com pu a 
tions were made from Eq. (3.3), (3.4) andTa S 
Estimates of actual concentrat ons may be dete ' 
mined by multiplying ordinate values by Q/ u 
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Figure 3-5A. xu/Q with distance for various heights of emission (H) and limits to vertical dispersion (L), A stability. 



Estimates 



11 



10 




DISTANCE, km 



Figure 3-5B. ,vu/Q with distance for various heights of emission (H) and limits to vertical dispersion (L), B stability. 
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Figure 3-5C. xu/Q with distance for various heights of emission (H) and limits to vertical dispersion (L), C stability. 
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Figure 3-5D. xu/Q with distance for various heights of emission (H) and limits to vertical dispersion (L), D stability. 
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Figure 3-5E. xii/Q with distance for various heights of emission (H) and limits to vertical dispersion (U, E stability. 
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igure 3-5F. xu/Q with distance for various heights of emission (H) and limits to vertical dispersion (L), F stabifity. 
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From Table A-l (Appendix 3) when exp 
"= 0-346, y/<r y = 1.46 



From Figure 3-2, for stability B and x = 600 m, <r y 
= 92. Therefore y = (1.46) (92) = 134 meters. 
This is the distance of the 10~ ;| isopleth from the 
x-axis at a downwind distance of 600 meters. 



This can also be determined from: 



2 In* 



x (x,y,0;H) 



(3.8) 



The position corresponding to the downwind dis- 
tance and off-axis distance can then be plotted. 
After a number of points have been plotted, the 
concentration isopleth may be drawn (see problems 
8 and 26). Figures 3-6 and 3-7 give ground-level 
isopleths of xu/Q for various stabilities for sources 
at H = and H = 100 meters. For example, to 
locate the 10~ s g m~ 3 isopleth resulting from a 
ground-level source of 20 g sec" 1 under B stability 
conditions with wind speed 2 m sec" 1 , one must 
first determine the corresponding value of x u /Q since 
this is the quantity graphed in Figure 3-6. x u /Q = 
10-" x 2/20 = 10-'. Therefore the x u/Q isopleth 
in Figure 3-6B having a value of 10"' m~ 3 corre- 
sponds to a x isopleth with a value of 10~ :) g m~ fl . 

AREAS WITHIN ISOPLETHS 

Figure 3-8 gives areas within isopleths of ground- 
level concentration in terms of x u/Q for a ground- 
level source for various stability categories (Gifford, 
1962; Hilsmeier and Gifford, 1962). For the exam- 
ple just given, the area of the 10~ a g m" 13 isopleth 
" u/Q isopleth) is about 5 x 10 4 meter 2 . 



(10~ J n 



CALCULATION OF MAXIMUM 
GROUND-LEVEL CONCENTRATIONS 

Figure 3-9 gives the distance to the point of 
maximum concentration, x limj; , and the relative maxi- 
mum concentration, x u/Q milx , as a function of 
effective height of emission and stability class 
(Martin, 1966). This figure was prepared from 
graphs of concentration versus distance, as in Fig- 
ure 3-5. The maximum concentration can be deter- 
mined by finding x u/Q as a function of effective 
emission height and stability and multiplying by 
Q/u. In using Figure 3-9, the user must keep in 
mind that the dispersion at higher levels may differ 
considerably from that determined by the o-/s and 
ff a 's used here. As noted, however, since a y gener- 
ally decreases with height and u increases with 



height, the product u a y a z will not change appreci- 
ably. The greater the effective height, the more 
likely it is that the stability may not be the same 
from the ground to this height. With the longer 
travel distances such as the points of maximum 
concentrations for stable conditions (Types E or 
F), the stability may change before the plume 
travels the entire distance. 

REVIEW OF ASSUMPTIONS 

The preceding has been based on these as- 
sumptions, which should be clearly understood: 

(i) Continuous emission from the source or 
emission times equal to or greater than travel times 
to the downwind position under consideration, so 
that diffusion in the direction of transport may be 
neglected. 

(ii) The material diffused is a stable gas or 
aerosol (less than 20 microns diameter) which re- 
mains suspended in the air over long periods of time. 

(iii) The equation of continuity: 



Q = 




X u cly dz 



(3.9) 



'In" denotes natural logarithms, i.e., to the base e, 



is fulfilled, i.e., none of the material emitted is re- 
moved from the plume as it moves downwind and 
there is complete reflection at the ground. 

( iv) The mean wind direction specifies the 
x-axis, and a mean wind speed representative of 
the diffusing layer is chosen. 

(v) Except where specifically mentioned, the 
plume constituents are distributed normally in both 
the cross-wind and vertical directions. 

(vi) The (r's given in Figures 3-2 and 3-3 repre- 
sent time periods of about 10 minutes. 
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Figure 3-8. Area within isopleths for a ground-level source (from Hilsmeier and Gifford). 
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Chapter 4 EFFECTIVE HEIGHT OF EMISSION 



GENERAL CONSIDERATIONS 

In most problems one must estimate the effec- 
tive stack height, H, at which the plume becomes 
essentially level. Rarely will this height correspond 
to the physical height of the stack, h. If the plume 
is caught in the turbulent wake of the stack or of 
buildings in the vicinity of the stack, the effluent 
will be mixed rapidly downward toward the ground 
(aerodynamic downwash). If the plume is emitted 
free of these turbulent zones, a number of emission 
factors and meteorological factors influence the rise 
of the plume. The emission factors are: velocity 
of the effluent at the top of the stack, V H ; tempera- 
ture of the effluent at the top of the stack, T H ; and 
diameter of the stack opening, d. The meteorolog- 
ical factors influencing plume rise are wind speed, 
u; temperature of the air, T,,; shear of the wind 
speed with height, du/dz; and atmospheric sta- 
bility. No theory on plume rise takes into account 
all of these variables; even if such a theory were 
available, measurements of all of the parameters 
would seldom be available. Most of the equations 
that have been formulated for computing the ef- 
fective height of emission are semi- empirical. For a 
recent review of equations for effective height of 
emission see Moses, Strom, and Carson (1964). 

Moses and Strom (1961), having compared ac- 
tual and calculated plume heights by means of six 
plume rise equations, report "There is no one for- 
mula which is outstanding in all respects." The 
formulas of Davidson-Bryant (1949), Holland 
(1953), Bosanquet-Carey-Halton (1950), and Bo- 
saiiquet (1957) all give generally satisfactory re- 
sults in the test situations. The experiments con- 
ducted by Moses and Strom involved plume rise 
from a stack of less than 0.5 meter diameter, stack 
gas exit velocities less than 15 m sec" 1 , and effluent 
temperature not more than 35C higher than that 
of the ambient air. 

The equation of Holland was developed with 
experimental data from larger sources than those 
of Moses and Strom (stack diameters from 1.7 to 
4,3 meters and stack temperatures from 82 to 
204 C); Holland's equation is used in the solution 
of the problems given in this workbook. This equa- 
tion frequently underestimates the effective height 
of emission; therefore its use often provides a slight 
"safety" factor. 

Holland's equation is; 



u 



where; 

AH = the rise of the plume above the stack, m 



V H = stack gas exit velocity, m sec l 
d = the inside stack diameter, m 
u = wind speed, m sec' 1 
p = atmospheric pressure, mb 
T H = stack gas temperature, K 
T a = air temperature, K 

and 2.68 x 10~- ! is a constant having units of mb" J 
m- 1 . 

Holland (1953) suggests that a value between 
1.1 and 1.2 times the AH from the equation should 
be used for unstable conditions; a value between 
0.8 and 0.9 times the AH from the equation should 
be used for stable conditions. 

Since the plume rise from a stack occurs over 
some distance downwind, Eq. (4.1) should not be 
applied within the first few hundred meters of the 
stack. 

EFFECTIVE HEIGHT OF EMISSION AND 
MAXIMUM CONCENTRATION 

If the effective heights of emission were the 
same under all atmospheric conditions, the highest 
ground-level concentrations from a given source 
would occur with the lightest winds. Generally, 
however, emission conditions are such that the ef- 
fective stack height is an inverse function of wind 
speed as indicated in Eq. (4.1). The maximum 
ground-level concentration occurs at some inter- 
mediate wind speed, at which a balance is reached 
between the dilution due to wind speed and the 
effect of height of emission. This critical wind speed 
will vary with stability. In order to determine the 
critical wind speed, the effective stack height as a 
function of wind speed should first be determined. 
The maximum concentration for each wind speed 
and stability can then be calculated from Figure 
3-9 as a function of effective height of emission 
and stability. When the maximum concentration 
as a function of wind speed is plotted on log-log 
graph paper, curves can be drawn for each stability 
class; the critical wind speed corresponds to the 
point of highest maximum concentration on the 
curve (see problem 14). 

ESTIMATES OF REQUIRED STACK HEIGHTS 

Estimates of the stack height required to pro- 
duce concentrations below a given value may be 
made through the use of Figure 3-9 by obtaining 
solutions for various wind speeds. Use of this figure 
considers maximum concentrations at any distance 
from the source. 

In some situations high concentrations upon the 
property of the emitter are of little concern, but 



Effective Height 



31 



maximum concentrations beyond the property line 
are of the utmost importance. For first approxima- 
tions it can be assumed that the maximum concen- 
tration occurs where \/Tcr B = H and that at this 
distance the IT'S are related to the maximum con- 
centration by: 



( r v a., ---= 



Q 



0.117 Q 



u e 



(4.2) 



Knowing the source strength, Q, and the concen- 
tration not to be exceeded xmnw one can determine 
the necessary a- y <r z for a given wind speed. Figure 
4-1 shows a-y a z as a function of distance for the 
various stability classes. The value of v y cr x and a 
design distance, x,i (the distance beyond which x is 
less than some pre-determined value) , will deter- 
mine a point on this graph yielding a stability class 
or point between classes. The CT K for this stability 
(or point between stabilities) can then be deter- 
mined from Figure 3-3. The required effective stack 
height for this wind speed can then be approxi- 
mated by H = \/2 <7 7 , (see problem 15). Since Eq. 
(4.2) is an approximation, the resulting height 
should be used with Eq. (3.3) to ensure that the 
maximum concentration is sufficiently low. If 
enough is known about the proposed source to 
allow use of an equation for effective height of 
emission, the relation between AH and u can be 
determined. The physical stack height required at 
the wind speed for which H was determined is H 
AH. The same procedure, starting with the deter- 
mination of a-y a z , must be used with other wind 
speeds to determine the maximum required physical 
stack height (see problem 16). 

EFFECT OF EVAPORATIVE COOLING 

When effluent gases are washed to absorb cer- 
tain constituents prior to emission, the gases are 
cooled and become saturated with water vapor. 
Upon release of the gases from the absorption tower, 
further cooling due to contact with cold surfaces 
of ductwork or stack is likely. This cooling causes 
condensation of water droplets in the gas stream. 
Upon release of the gases from the stack, the water 
droplets evaporate, withdrawing the latent heat of 
vaporization from the air and cooling" the plume. 
The resulting negative buoyancy reduces the effec- 
tive stack height (Scorer, 1959). 

EFFECT OF AERODYNAMIC DOWNWASH 

The influence of mechanical turbulence around 
a building or stack can significantly alter the ef- 
fective stack height. This is especially true with 
high winds, when the beneficial effect of high stack- 
gas velocity is at a minimum and the plume is 
emitted nearly horizontally. The region of disturbed 
flow surrounds an isolated building, generally to at 



least twice its height and extends downwind 5 to 10 
times its height. Building the stack 2.5 times the 
height of the highest building adjacent to the stack 
usually overcomes the effects of building turbulence 
(Hawkins and Nonhebel, 1955). Ensuring that the 
exit velocity of the stack gas is more than 1.5 times 
the wind speed will usually prevent downwash in 
the wake of the stack. Most of the knowledge about 
the turbulent wakes around stacks and buildings 
has been gained through wind tunnel studies (Sher- 
lock and Lesher, 1954; Strom, 1955-1956; Strom, 
et al, 1957; and Halitsky, 1962). By use of models 
of building shapes and stacks, one may determine 
the wind speeds required to cause downwash for 
various wind directions. With a wind tunnel the 
meteorological variables most easily accounted for 
are wind speed and wind direction (by rotation of 
the model within the tunnel) . The emission factors 
that may be considered are the size and shape of 
the plant building; the shape, height, and diameter 
of the stack; the amount of emission; and the stack- 
gas velocity. 

Through wind tunnel studies, the critical wind 
speeds that will cause downwash from various di- 
rections can be determined for a given set of plant 
factors. The average number of hours of downwash 
per year can then be calculated' by determining the 
frequency of wind speeds greater than the critical 
speeds for each direction (Sherlock and Lesher, 
1954) if climatological data representative of the 
site are available. 

Maximum downwash about a rectangular struc- 
ture occurs when the direction of the wind is at an 
angle of 45 degrees from the major axis of the struc- 
ture; minimum downwash occurs with wind flow 
parallel to the major axis of the structure (Sherlock 
and Lesher, 1954). 

Halitsky (1961, 1963) has shown that the efflu- 
ent from flush openings on flat roofs frequently 
flows in a direction opposite to that of the free 
atmospheric wind, owing to counter-flow along the 
roof in the turbulent wake above the building. In 
addition to the effect of aerodynamic downwash 
upon the release of air pollutants from stacks and 
buildings, one must also consider the effects of aero- 
dynamic downwash when exposing meteorological 
instruments near or upon buildings. 

Where the pollution is emitted from a vent or 
opening on a building and is immediately influ- 
enced by the turbulent wake of the building, the 
pollution is rapidly distributed within this turbu- 
lent wake. To account for mixing in the turbulent 
wake, one may assume binomial distributions of 
concentrations at the source, with horizontal and 
vertical standard deviations of <r yo and a m . The 
standard deviations are related to the width and 
height of the building, for example, letting 4.3 u yf > 
equal the width of the building and 2.15 a. M equal 
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Figure 4-1. The product of <vr. as a function of downwind distance from the source. 
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the height. Values other than 4.3 and 2.15 can be 
used. When these values are used 97 % of the dis- 
tribution is included within these limits. Virtual 
distances x y and X K can be found such that at x -v , 
ff y IT.,.,, and at x a , a x , == fr 7n . These x's will differ 
with stability. Equations applicable to point sources 
can then he used, determining <r v as a function of 
x -f~ x,. and T Z as a function of x -f x ?p . 
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CONCENTRATIONS IN AN INVERSION 
BREAK-UP FUMIGATION 

A surface-based inversion may be eliminated by 
the upward transfer of sensible heat from the 
ground surface when that surface is warmer than 
the overlying air. This situation occurs when the 
ground is being warmed by solar radiation or when 
air flows from a cold to a relatively warm surface. 
In either situation pollutants previously emitted 
above the surface into the stable layer will be mixed 
vertically when they are reached by the thermal 
eddies, and ground-level concentrations can increase. 
This process, called "fumigation" was described by 
Hewson and Gill (1944) and Hewson (1945). Equa- 
tions for estimating concentrations with these con- 
ditions have been given by Holland (1953), Hew- 
son (1955), Gilford (1960a), Bierly and Hewson 
(1962), and Pooler (1965). 

To estimate ground-level concentrations under 
inversion break-up fumigations, one assumes that 
the plume was initially emitted into a stable layer. 
Therefore, <r y and cr K characteristic of stable ^ condi- 
tions must be selected for the particular distance 
of concern. An equation for the ground-level con- 
centration when the inversion has been eliminated 
to a height hi is: 



(x,y,0;H) = 



p 



Q 



exp ( 0.5 p a ) dp _ 



exp 



STT ffyif U hf 

i y 



(5.1) 



h, H 



where p 

"z 

and (TyF is discussed below. 

Values for the integral in brackets can be found in 
most statistical tables. For example, see pages 273- 
276, Burington (1963). This factor accounts for 
the portion of the plume that is mixed downward. 
If the inversion is eliminated up to the effective 
stack height, half of the plume is presumed to be 
mixed downward, the other half remaining in the 
stable air above. Eq. (5.1) can be approximated 
when the fumigation concentration is near its 
maximum by: 



(x,y,0;H) 




Q 



(5.2) 
(5.3) 



A difficulty is encountered in estimating a rea- 
sonable value for the horizontal dispersion since in 
mixing the stable plume through a vertical depth 
some additional horizontal spreading occurs (see 
problem 12). If this spreading is ignored and the 
o-j- for stable conditions used, the probable result 
would be estimated concentrations higher than ac- 
tual concentrations. Or, using an approximation 
suggested by Bierly and Hewson (1962) that the 
edge of the plume spreads outward with an angle 
of 15, the cfyit for the inversion break-up fumigation 
equals the a y for stable conditions plus one-eighth 
the effective height of emission. The origin of this 
concept can be seen in Figure 5-1 and the following 
equation, where the edge of the plume is the point 
at which the concentration falls to 1/10 that at the 
centerline (at a distance of 2.15 * y from the plume 
center) , 

2.15 ffs . (stable) + H tan 15 



" rff 2.15 

-= <r y (stable) -f H/8 (5.4) 

A Gaussian distribution in the horizontal is as- 
sumed. 
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Figure 5-1. Diagram showing assumed height, hi and *? 
during fumigation, for use in equation (5.2). 

Eq. (6.4) should not be applied near the stack, 
for if the inversion has been eliminated to a height 
sufficient to include the entire plume, the emission 
is taking place under unstable not stable conditions. 
Therefore, the nearest downwind distance to be 
considered for an estimate of fumigation concen- 
trations must be great enough, based on the time 
required to eliminate the inversion, that this por- 
tion of the plume was initially emitted into stable 
air. This distance is x =- ut m , where u is the mean 
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wind in the stable layer and t rii is the time required 
to eliminate the inversion from h, the physical 
height of the stack to h f (Eq. 5.3). 

t,,, is dependent upon both the strength of the 
inversion and the rate of heating at the surface. 
Pooler (1965) has derived an expression for esti- 
mating this time: 



p:i 



SO 



R 



(hi h) 



h, 



(5.5) 



where t m 



time required for the mixing layer to 
develop from the top of the stack to the 
top of the plume, sec 
p a = ambient air density, g m~ 3 

Cp = specific heat of air at constant pressure, 
cal g- 1 "K- 1 

R =-- net rate of sensible heating of an air 
column by solar radiation, cal rn~ 2 sec" 1 

so 

vertical potential temperature gradient, 



"Km- 1 
rate) 



Sz 



h r (the adiabatic lapse 



hi = height of base of the inversion sufficient 
to be above the plume, m 

h = physical height of the stack, m 
Note that hi h is the thickness of the layer to be 
heated and f - 1 1 is the average height of the 

layer. Although R depends on season, and cloud 
cover and varies continuously with time, Pooler has 
used a value of 67 cal m~ 2 sec" 1 as an average for 
fumigation. 

Hewson (1945) also suggested a method of esti- 
mating the time required to eliminate an inversion 
to a height z by use of an equation of Taylor's 
(1915, p. 8): 



u ~ 4 K 

where: t = time required to eliminate the inver- 
sion to height z, sec 

z = height to which the inversion has been 
eliminated, m 

K = eddy diffusivity for heat, m 2 sec" 1 
Rewriting to compare with Eq. (5.5), 
V h 2 



4 K 



(5.7) 



Hewson (1945) has suggested a value of 3 m 3 sec" 1 
for K. 

PLUME TRAPPING 

Plume trapping occurs when the plume is 
trapped between the ground surface and a stable 



layer aloft. Bierly and Hewson (1962) have sug- 
gested the use of an equation that accounts for the 
multiple eddy reflections from both the ground and 
the stable layer: 
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where L is the height of the stable layer and J = 3 
or 4 is sufficient to include the important reflec- 
tions. A good approximation of this lengthy equa- 
tion can be made by assuming no effect of the stable 
layer until <r a = 0.47 L (see Chapter 3). It is as- 
sumed that at this distance, XL, the stable layer 
begins to affect the vertical distribution so that at 
the downwind distance, 2 XL, uniform vertical mix- 
ing has taken place and the following equation can 
be used: 
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For distances between XL and 2 X L the best approxi- 
mation to the ground-level centerline concentration 
is that read from a straight line drawn between the 
concentrations for points XL and 2 XL on a log-log 
plot of ground-level centerline concentration as a 
function of distance. 

CONCENTRATIONS AT GROUND LEVEL 
COMPARED TO CONCENTRATIONS AT THE 
LEVEL OF EFFECTIVE STACK HEIGHT 
FROM ELEVATED CONTINUOUS SOURCES 

There are several interesting relationships be- 
tween ground-level concentrations and concentra- 
tions at the level of the plume centerline. One of 
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these is at the distance of maximum concentration 
at the ground. As a rough approximation the maxi- 
mum ground-level concentration occurs at the dis- 
1 



tance where 



H. This approximation is 



much better for unstable conditions than for stable 
conditions. With this approximation, the ratio of 
concentration at plume centerline to that at the 
ground is: 
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This calculation indicates that at the distance 
of maximum ground-level concentration the concen- 
tration at plume centerline is greater by about 
one- third. 

It is also of interest to determine the relation- 
ship between <r K and H such that the concentration 
at ground-level at a given distance from the source 
is the same as the concentration at plume level. 
This condition should occur where: 



1 / H 



1.0 + exp 



The value HA K = 1.10 satisfies this expression, 
which can be written as <r z = 0.91 H (see problem 
10). 

TOTAL DOSAGE FROM A FINITE RELEASE 

The total dosage, which is the integration of 
concentration over the time of passage of a plume 
or puff, can be obtained from; 



D* (x,y,0;H) = 
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where D v = total dosage, g sec m 3 
and Q T = total release, g 



The ff's should be representative of the time 
period over which the release takes place, and care 
should be taken to consider the x-axis along the 
trajectory or path of the plume or puff travel. Large 
errors can easily occur if the path is not known 



accurately. The estimate of this path is usually in- 
creasingly difficult with shorter release times. D T 
can also be given in curie sec m~ 3 if QT is in curies. 

CROSSWIND-INTEGRATED CONCENTRATION 

The ground-level cross wind-in teg rated concen- 
tration is often of interest. For a continuous ele- 
vated source this concentration is determined from 
Eq. (3.2) integrated with respect to y from ~ x to 
+ " (Gifford 1960a) giving: 
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In diffusion experiments the ground-level cross- 
wind-integrated concentration is often determined 
at particular downwind distances from a crosswind 
line or arc of sampling measurements made at this 
distance. When the source strength, Q, and average 
wind speed, u, are known, o x can be estimated in- 
directly even though no measurements were made 
in the vertical. If any of the tracer is lost through 
reaction or deposition, the resulting a z from such 
estimates will not represent the vertical dispersion 
(see problem IS). 

ESTIMATION OF CONCENTRATIONS FOR 
SAMPLING TIMES LONGER THAN A 

FEW MINUTES 

Concentrations directly downwind from a source 
decrease with sampling time mainly because of a 
larger <r y due to increased meander of wind direction. 
Stewart, Gale, and Crooks (1958) reported that 
this decrease in concentration follows a one-fifth 
power law with the sampling time for sampling 
periods from about 3 minutes to about half an hour. 
Cramer (1959) indicates that this same power law 
applies for sampling times from 3 seconds to 10 
minutes. Both of these studies were based on ob- 
servations taken near the height of release. Gifford 
(1960b) indicates that ratios of peak to mean con- 
centrations are much higher than those given by 
the above power law where observations of concen- 
trations are made at heights considerably different 
from the height of release or considerably removed 
from the plume axis. He also indicates that for 
increasing distances from an elevated source, the 
ratios of peak to average concentrations observed 
at ground level approach unity. Singer (1961) and 
Singer, et al. (1963) show that ratios of peak to 
mean concentrations depend also- on the stability 
of the atmosphere and the type of terrain that the 
plume is passing over. Nonhebel (I960) reports 
that Meade deduced a relation between calculated 
concentrations at ground level and the sampling 
time from "a study of published data on lateral and 
vertical diffusion coefficients in steady winds." 
These relations are shown in Table 5-1, 



Special Topics 



37 



Table 5-1 VARIATION OF CALCULATED CONCENTRATION 
WITH SAMPLING TIME 

Ratio of 

Calculated Concentration 



2 Q 



Sampling Time 


to 3-minute Concentration 


3 minutes 


1.00 


15 minutes 


0.82 


1 hour 


0.61 


3 hours 


0.51 


24 hours ... 


0.36 



This table indicates a power relation with time: 
X <x t~ ft - 1T . Note that these estimates were based 
upon published dispersion coefficients rather than 
upon sampling results. Information in the refer- 
ences cited indicates that effects of sampling time 
are exceedingly complex. If it is necessary to esti- 
mate concentrations from a single source for the 
time intervals greater than a few minutes, the best 
estimate apparently can be obtained from; 



X = 



(5.12) 



vhere ,\ is the desired concentration estimate for 
;he sampling time, t s ; x* is the concentration esti- 
nate for the shorter sampling time, t k> (probably 
ibout 10 minutes); and p should be between 0.17 
md 0.2. Eq. (5.12) probably would be applied 
nost appropriately to sampling times less than 2 
lours (see problem 19). 

ESTIMATION OF SEASONAL OR ANNUAL 
AVERAGE CONCENTRATIONS AT A 
RECEPTOR FROM A SINGLE POLLUTANT 
SOURCE 

For a source that emits at a constant rate from 
lour to hour and day to day, estimates of seasonal 
jr annual average concentrations can be made for 
any distance in any direction if stability wind "rose" 
data are available for the period under study. A 
wind rose gives the frequency of occurrence for 
each wind direction (usually to 16 points) and wind 
speed class (9 classes in standard Weather Bureau 
use) for the period under consideration (from 1 
month to 10 years) . A stability wind rose gives the 
same type of information for each stability class. 

If the wind directions are taken to 16 points and 
it is assumed that the wind directions within each 
sector are distributed randomly over a period of a 
month or a season, it can further be assumed that 
the effluent is uniformly distributed in the hori- 
zontal within the sector (Holland, 1953, p. 540). 
The appropriate equation for average concentration 
is then either: 
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or 
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(5.14) 



depending upon whether a stable layer aloft is af- 
fecting the distribution. 

The estimation of x for a particular direction 
and downwind distance can be accomplished by 
choosing a representative wind speed for each speed 
class and solving the appropriate equation (5.13 or 
5.14) for all wind speed classes and stabilities. Noto 
that a SSW wind affects a receptor to the NNE 
of a source. One obtains the average concentration 
for a given direction and distance by summing all 
the concentrations and weighting each one accord- 
ing to its frequency for the particular stability and 
wind speed class. If desired, a different effectivo 
height of emission can be used for various wind 
speeds. The average concentration can be expressed 
by: 
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where f (0, S, N) is the frequency during the period 
of interest that the wind is from the direc- 
tion G, for the stability condition, S, and 
wind speed class N. 

o- zS is the vertical dispersion parameter evaluated 
at the distance x for the stability condition S. 

UN is the representative wind speed for class N, 

H,, is the effective height of release for the wind 
speed UN. 

Where stability wind rose information cannot bo 
obtained, a first-order approximation may be made 
of seasonal or annual average concentrations by 
using the appropriate wind rose in the same man- 
ner, and assuming the neutral stability class, D, 
only. 

METEOROLOGICAL CONDITIONS 
ASSOCIATED WITH MAXIMUM 
GROUND-LEVEL CONCENTRATIONS 

1, For ground-level sources maximum concentra- 
tions occur with stable conditions. 
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2. For elevated sources maximum "instantaneous" 
concentrations occur with unstable conditions 
when portions of the plume that have undergone 
little dispersion are brought to the ground. 
These occur close to the point of emission (on 
the order of 1 to 3 stack heights) . These con- 
centrations are usually of little general interest 
because of their very short duration; they can- 
not be estimated from the material presented in 
this workbook. 

3. For elevated sources maximum concentrations 
for time periods of a few minutes occur with 
unstable conditions; although the concentra- 
tions fluctuate considerably under these condi- 
tions, the concentrations averaged over a few 
minutes are still high compared to those found 
under other conditions. The distance of this 
maximum concentration occurs near the stack 
(from 1 to 5 stack heights downwind) and the 
concentration drops off rapidly downwind with 
increasing distance, 

4. For elevated sources maximum concentrations 
for time periods of about half an hour can occur 
with fumigation conditions when an unstable 
layer increases vertically to mix downward a 
plume previously discharged within a stable 
layer. With small AH, the fumigation can occur 
close to the source but will be of relatively short 
duration. For large AH, the fumigation will 
occur some distance from the stack (perhaps 30 
to 40 km), but can persist for a longer time 
interval. Concentrations considerably lower than 
those associated with fumigations, but of sig- 
nificance can occur with neutral or unstable 
conditions when the dispersion upward is se- 
verely limited by the existence of a more stable 
layer above the plume, for example, an inversion. 

5. Under stable conditions the maximum concen- 
trations at ground-level from elevated sources 
are less than those occurring under unstable 
conditions and occur at greater distances from 
the source. However, the difference between 
maximum ground-level concentrations for stable 
and unstable conditions is only a factor of 2 
for effective heights of 25 meters and a factor 
of 6 for H of 75 m. Because the maximum 
occurs at greater distances, concentrations that 
are below the maximum but still significant can 
occur over large areas. This becomes increas- 
ingly significant if emissions are coming from 
more than one source. 

CONCENTRATIONS AT A RECEPTOR POINT 
FROM SEVERAL SOURCES 

Sometimes, especially for multiple sources, it is 
convenient to consider the receptor as being at the 
origin of the diffusion coordinate system. The 



source-receptor geometry can then be worked out 
merely by drawing or visualizing an x-axis oriented 
upwind from the receptor and determining the 
crosswind distances of each source in relation to this 
x-axis. As pointed out by GifEord (1959), the con- 
centration at (0, 0, 0) from a source at (x, y, H) 
on a coordinate system with the x-axis oriented up- 
wind is the same as the concentration at (x, y, 0) 
from a source at (0, 0, H) on a coordinate system 
with the x-axis downwind (Figure 5-2), The total 
concentration is then given by summing the indi- 
vidual contributions from each source {see problem 
20). 

SOURCE 

UPWIND 



RECEPTOR 
(0,0,0) 




DOWNWEND 



Figure 5-2. Comparison of source-oriented and receptor- 
oriented coordinate systems. 

It is often difficult to determine the atmos- 
pheric conditions of wind direction, wind speed, and 
stability that will result in the maximum combined 
concentrations from two or more sources; drawing 
isopleths of concentration for various wind speeds 
and stabilities and orienting these according to 
wind direction is one approach. 

AREA SOURCES 

In dealing with diffusion of air pollutants in 
areas having large numbers of sources, e.g., as in 
urban areas, there may be too many sources of most 
atmospheric contaminants to consider each source 
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individually. Often an approximation can be made 
by combining all of the emissions in a given area 
and treating this area as a source having an initial 
horizontal standard deviation, tr voi A virtual dis- 
tance, x y , can then be found that will give this 
standard deviation. This is just the distance that 
will yield the appropriate value for a y from Figure 
3-2. Values of x, will vary with stability. Then 
equations for point sources may be used, determin- 
ing ff y as a function of x -{- x y , a slight variation of 
the suggestion by Holland {1953). This procedure 
treats the area source as a cross-wind line source 
with a normal distribution, a fairly good approxi- 
mation for the distribution across an area source. 
The initial standard deviation for a square area 
source can be approximated by <r yo , s/4.3, where 
s is the length of a Bide of the area (see problem 
22). 

If the emissions within an area are from varying 
effective stack heights, the variation may be ap- 
proximated by using a a. Thus H would be the 
mean effective height of release and o zo the standard 
deviation of the initial vertical distribution of 
sources, A virtual distance, x z , can be found, and 
point source equations used for estimating concen- 
trations, determining <j y . as a function of x + x z . 

TOPOGRAPHY 

Under conditions of irregular topography the 
direct application of a standard dispersion equation 
is often invalid. In some situations the best one 
may be able to do without the benefit of in situ 
experiments is to estimate the upper limit of the 
concentrations likely to occur. 

For example, to calculate concentrations on a 
hillside downwind from and facing the source and 
at about the effective source height, the equation 
for concentrations at ground-level from a ground- 
level source (Eq* 3,4) will yield the highest ex- 
pected concentrations. This would closely approxi- 
mate the situation under stable conditions, when 
the pollutant plume would be most likely to en- 
counter the hillside. Under unstable conditions the 
flow is more likely to rise over the hill (see problem 
21). 

With downslope flow when the receptor is at a 
lower elevation than the source, a likely assumption 
is that the flow parallels the slope; i.e., no allow- 
ance is made for the difference between ground- 
level elevations at the source and at the receptor. 

Where a steep ridge or bluff restricts the hori- 
zontal dispersion, the flow is likely to be parallel 
to such a bluff. An assumption of complete reflec- 
tion at the bluff, similar to eddy reflection at the 
ground from an elevated source, is in order. This 
may be accomplished by using: 
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B is the distance from the x-axis to the restrict- 
ing bluff, and the positive y axis is defined to be in 
the direction of the bluff. 

The restriction of horizontal dispersion by valley 
sides is somewhat analogous to restriction of the 
vertical dispersion by a stable layer aloft. When 
the fr y becomes great enough, the concentrations 
can be assumed to be uniform across the width of 
the valley and the concentration calculated accord- 
ing to the following equation, where in this case Y 
is the width of the valley. 
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LINE SOURCES 



Concentrations downwind of a continuously 
emitting infinite line source, when the wind direc- 
tion is normal to the line, can be expressed by 
rewriting equation (12) p. 154 of Sutton (1932): 

X (x,y,0;H) = . q exp 1 ^ 

VSwr cr a U L 

(5.18) 

Here q is the source strength per unit distance, 
for example, g sec" 1 m -1 . Note that the horizontal 
dispersion parameter, a y , does not appear in this 
equation, since it is assumed that lateral dispersion 
from one segment of the line is compensated by dis- 
persion in the opposite direction from adjacent 
segments. Also y does not appear, since concentra- 
tion at a given x is the same for any value of y 
(see problem 23). 

Concentrations from infinite line sources when 
the wind is not perpendicular to the line can be 
approximated. If the angle between the wind direc- 
tion and line source is 0, the equation for concen- 
tration downwind of the line source is: 



X (x,y,0;H) 
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This equation should not be used where is les* 
than 45. 
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When estimating concentrations from finite line 
sources, one must account for "edge effects" caused 
by the end of the line source. These effects will of 
course extend to greater cross-wind distances as 
the distance from the source increases. For concen- 
trations from a finite line source oriented cross- 
wind, define the x-axis in the direction of the mean 
wind and passing through the receptor of interest. 
The limits of the line source can be defined as ex- 
tending from y, to y, where y, is less than y a . The 
equation for concentration (from Button's (1932) 
equation (11), p. 154), is: 



(x,0,0;H) 




(5.20) 



The value of the integral can be determined from 
tabulations given in most statistical tables (for ex- 
ample, see Burrington (1953), pp. 273-276; also see 
problem 24). 

INSTANTANEOUS SOURCES 

Thus far we have considered only sources that 
were emitting continuously or for time periods equal 
to or greater than the travel times from the source 
to the point of interest. Cases of instantaneous re- 
lease, as from an explosion, or short-term releases 
on the order of seconds, are often of practical con- 
cern. To determine concentrations at any position 
downwind, one must consider the time interval 
after the time of release and diffusion in the down- 
wind direction as well as lateral and vertical diffu- 
sion. Of considerable importance, but very difficult, 
is the determination of the path or trajectory ot 
the "puff." This is most important if concentra- 
tions are to be determined at specific points. Deter- 
mining tbe trajectory is of less importance if knowl- 
edge of the magnitude of the concentrations tor 
particular downwind distances or travel times is 
required without the need to know exactly at what 
points these concentrations occur. Rewriting tout- 
ton's (1932) equation (13), p. 155, results in an 
equation that may be used for estimates of concen- 
tration downwind from a release from height, n. 
. . 2 QT r 1 
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The symbols have the usual meaning, with the 
important exceptions that Q T represents the total 
mass of the release and the u's are not those eval- 
uated with respect to the dispersion of a continuous 
source at a fixed point in space. 

In Eq. (5.21) the </s refer to dispersion sta- 
tistics following the motion of the expanding puff. 
The v x is the standard deviation of the concentra- 
tion distribution in the puff in the downwind direc- 
tion, and t is the time after release. Note that 
there is no dilution in the downwind direction by 
wind speed. The speed of the wind mainly serves 
to give the downwind position of the center of the 
puff, as shown by examination of the exponential 
involving cr s . Wind speed may influence the dis- 
persion indirectly because the dispersion parameters 
CT XJ try, and <J B may be functions of wind speed. The 
^'s and ff,'s for an instantaneous source are less 
than those for a few minutes given in Figure 3-2 and 
3-3. Slade (1965) has suggested values for a a, 
and a z for quasi-instantaneous sources. These are 
given in Table 5-2. The problem remains to make 
best estimates of a x . Much less is known of diffu- 
sion in the downwind direction than is known oE 
lateral and vertical dispersion. In general one should 
expect the ^ value to be about the same as o f . 
Initial dimensions of the puff, i.e., from an explo- 
sion, may be approximated by finding a virtual 
distance to give the appropriate initial standard 
deviation for each direction. Then * y will be deter- 
mined as a function of x + * * as a function of 
x + X H and ^ as a function of x + x*. 

Table 5-2 ESTIMATION OF DISPERSION PARAMETERS FOR 
QUASI-INSTANTANEOUS SOURCES (FROM SLADE. 1965) 




Unstable 
Neutral 
Very Stable 
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Chapter 6 RELATION TO OTHER DIFFUSION EQUATIONS 



Most other widely used diffusion equations are 
variant forms of the ones presented here. With re- 
spect to ground-level concentrations from an ele- 
vated source (Eq. 3.2): 



(x,y,0;H) = 
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Other well-known equations can be compared: 
Bosanquet and Pearson (1936): 
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where p and q are dimensionless diffusion coeffi- 
cients. 
Sutton (1947): 
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where n is a dimensionless constant and C y and C a 
are diffusion coefficients in m"' 2 . 

Calder (1952): 
x <x, yj O;H) = 



Q u 



2 k 8 a v, 2 x 2 



exp 



4- H 



(6.3) 



where a = -^7 , the ratio of horizontal eddy velocity 

w' 

to vertical eddy velocity, k is von Karman's con- 
stant approximately equal to 0.4, and v x = g 

, A In ( ) 

where z u is a roughness parameter, m. z 



NOTE: Calder wrote the equation for the con- 
centration at (x, y, z) from a ground-level source. 
For Eq. (6.3) it is assumed that the concentration 
at ground level from an elevated source is the same 
as the concentraton at an elevated point from a 
ground-level source. 

Table 6-1 lists the expressions used in these 
equations that are equivalent to <r y and <T K (con- 
tinuous source) in this paper. 

Table 6-1 EXPRESSIONS EQUIVALENT TO ^ AND , IN 
VARIOUS DIFFUSION EQUATIONS. 
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Chapter 7 EXAMPLE PROBLEMS 



The following 26 example problems and their 
solutions illustrate the application of most of the 
techniques and equations presented in this work- 
book. 

PROBLEM 1: It is estimated that a burning 
dump emits 3 g sec" 1 of oxides of nitrogen. 
What is the concentration of oxides of nitrogen, 
averaged over approximately 10 minutes, from 
this source directly downwind at a distance of 
3 km on an overcast night with wind speed ^ of 
7 m sec" 1 ? Assume this dump to be a point 
ground-level source with no effective rise. 

SOLUTION: Overcast conditions with a wind 
speed of 7 m sec" 1 indicate that stability class D 
is most applicable (Statement, bottom of Table 
3-1). For x 3 km and stability D, o> = 190 m 
from Figure 3-2 and <r a = 65 m from Figure 3-3. 
Eq. (3.4) for estimation of concentrations di- 
rectly downwind (y = 0) from a ground-level 
source is applicable: 



X <x,0,G;0) - 



Q 



U 



TT 190 (65) 7 



= 1.1 x 10"* g m~ :t of oxides of nitrogen. 

PROBLEM 2: It is estimated that 80 g sec" 1 of 
sulfur dioxide is being emitted from a petroleum 
refinery from an average effective height of 60 
meters. At 0800 on an overcast winter morning 
with the surface wind 6 m sec" 1 , what is the 
ground-level concentration directly downwind 
from the refinery at a distance of 500 meters? 

SOLUTION: For overcast conditions, D class sta- 
bility applies. With D stability at x = 600 m, 
ffy 36 m, ffi = 18.5 m. Using Eq. (3.3): 

Q r i / H " 



X (x,0,0;H) = 



U 



exp 



80 



7r36 (18.5) < 
-= 6.37 x 10" 3 exp [0.5 (3.24) 3 ] 

The exponential is solved using Table A-l (Ap- 
pendix 3). 

==6.37xlO" 3 (5.25 xlO" 3 ) 
X 3.3 x 10~ 6 g m" a of S0 2 

PROBLEM 3: Under the conditions of problem 
2, what is the concentration at the same dis- 
tance downwind but at a distance 50 meters 
from the x-axis? That is: x (500, 50, 0; 60) =- ? 

SOLUTION: Using Eq. (3.2): 

Q ..... f i ( y 



(x,y,0;H) 



U 



exp 



exp 



H 



All but the exponential involving y has been 
found in the preceding problem. Therefore: 

x (500, 50,0; 60) =3.3x10-' 
exp [0.5 (50/36)*] 

= 3.3 xlQ~ B (0,381) 
'= 1.3 x 10- fi g m~ 3 of S0 a 

PROBLEM 4: A power plant burns 10 tons per 
hour of coal containing 3 percent sulfur; the 
effluent is released from a single stack. On a 
sunny summer afternoon the wind at 10 meters 
above ground is 4 m sec" 1 from the northeast. 
The morning radiosonde taken at a nearby 
Weather Bureau station has indicated that a 
frontal inversion aloft will limit the vertical 
mixing to 1500 meters, The 1200-meter wind is 
from 30 at 5 m sec" 1 . The effective height of 
emission is 150 meters. From Figure 3-9, what 
is the distance to the maximum ground-level 
concentration and what is the concentration at 
this point? 

SOLUTION: To determine the source strength, 
the amount of sulfur burned is; 10 tons hr" 1 x 
2000 Ib ton-' x 0.03 sulfur = 600 Ib sulfur hr 1 . 
Sulfur has a molecular weight of 32 and com- 
bines with 3 with a molecular weight of 32; 
therefore for every mass unit of sulfur burned, 
there result two mass units of S0 2 . 



Q 



X 



64 (molecular weight of S0 a ) 
32 (molecular weight of sulfur) 
eOQlbhr" 1 (453.6 gib" 1 ) 

3600 sec hr" 1 
ISlgsec-'ofSO, 



On a sunny summer afternoon the insolation 
should be strong. From Table 3-1, strong inso- 
lation and 4m sec" 1 winds yield class-B stability. 
From Figure 3-9, the distance to the point of 
maximum concentration is 1 km for class-B sta- 
bility and effective height of 150 meters. From 
Figure 3-3 at this distance <r a = 110 m. This is 
much less than 0.47 L. Therefore, at this dis- 
tance, the limit of mixing of 1500 meters will 
not affect the ground-level concentration. From 
Figure 3-9, the maximum *u/Q for B stability 
and this effective height of 150 m is 7.5 x 10 . 



Q 



7.5 x 10" x 151 



u 



= 2.8 x 10"* g m" 3 of S0 a 

PROBLEM 5: For the power plant in problem 4, 
at what distance does the maximum ground- 
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level concentration occur and what is this con- 
centration on an overcast day with wind speed 
4 m sec" 1 ? 

SOLUTION: On an overcast day the stability 
class would be D. From Figure 3-9 for D sta- 
bility and H of 150 m, the distance to the point 
of maximum ground-level concentration is 5.6 
km, and the maximum xu/Q is 3.0 x 1Q~ U . 

3.0 x lET* x 151 



= 1.1 x ID- 1 g m" 3 

PROBLEM 6: For the conditions given in prob- 
lem 4, draw a graph of ground-level centerline 
sulfur dioxide concentration with distance from 
100 meters to 100 km. Use log-log graph paper. 

SOLUTION: The frontal inversion limits the mix- 
ing to L = 1500 meters. The distance at which 
a, = Q.47 L = 705 m is XL = 5.6 km. At dis- 
tances less than this, Eq. (3.3) is used to calcu- 
late concentrations: 



x (x,0,0;H) = 



Q 



exp 



At distance equal to or greater than 2 XL, which 
is 11 km, Eq, (3.5) is used: 



(x,0,0;H) = 



Q 



2?r CTV L U 



Solutions for the equations are given in Table 
7-1. The values of concentration are plotted 
against distance in Figure 7-1. 




DOWNWIND DISTANCI, km 



Figure 7-1. Concentration as a lunction of downwind 
distance (Problem 6). 



Table 7-1 CALCULATION OF CONCENTRATIONS FOR 
VARIOUS DISTANCES (PROBLEM 6) 



x, 
km 


u, 
m sec" 1 


<fy, 

m 


m 


H/*. 


p[-f<H/ ffl ; 


"1. 


X- 
nr Ll 


0. 


3 


4 




52 


30 


5.0 


3.73 x 10~ 


2.9 


X 


10- 


s 


0. 


5 


4 




83 


51 


2.94 


1.33 x 10-" 


3.8 


X 


10- 


a 


0.8 


4 




129 


85 


1.77 


0.209 


2.3 


X 


10- 


1 


1. 





4 




157 


no 


1.36 


0.397 


2.8 


X 


10 


1 


2. 





4 




295 


230 


0.65 


0.810 


1.4 


X 


10- 


4 


3. 





4 




425 


365 


0.41 


0.919 


7.1 


X 


10- 


a 


5. 


5 


4 


.5 


720 


705 


0.21 


0.978 


2.1 


X 


10" 


.1. 


x, 
km 


Hi 

m sec" 1 


ffy, 

m 


L, 
m 






B 


X- 
Ml ' 


11. 





4 


.5 


1300 


1500 






6.9 


X 


10 


a 


30 


4.5 


3000 


1500 






3.0 


X 


10 


o 


100 


4.5 


8200 


1500 






1.1 


X 


10- 






PROBLEM 7: For the conditions given in prob- 
lem 4, draw a graph of ground-level concentra- 
tion versus crosswind distance at a downwind 
distance of 1 km. 

SOLUTION: From problem 4 the ground-lovcl 
centerline concentration at 1 km is 2.8 x 10'"* 
g m~ 3 . To determine the concentrations at dis- 
tances y from the x-axis, the ground-level con- 
terline concentration must be multiplied by tho 



factor exp 



o-y = 157 meters at x = 1 km. Values for this 
computation are given in Table 7-2. 

Table 7-2 DETERMINATION OF CROSSWIND 
CONCENTRATIONS (PROBLEM 7) 



y, 

m 


y 


exo f 1 f-! 

exp l *U 


)1 


x (x,y r O) 


100 


0.64 


0.815 




2.3x10'' 


200 


1.27 


0.446 




1.3 xlO" 1 


300 


1.91 


0.161 




4r i /i -^ FI 
.5x 10 a 


400 


2.55 


3.87 x 10- 3 




i.ixio- s 


500 


3.18 


6.37 x 1Q- 3 




1.8x10 " 



These concentrations are plotted in Figure 7-2. 

PROBLEM 8: For the conditions given in prob- 
lem 4, determine the position of the 10" fl m"* 
ground level isopleth, and determine its urea, 

SOLUTION: From the solution to problem 6, the 
graph (Figure 7-1) shows that the 1Q- B g irr 3 
isopleth intersects the x-axis at approximately 
x = 350 meters and x =- 8.6 kilometers. 
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-200 +200 

CROSSWIND DISTANCE ly), m 
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Figure 7-2. Concentration as a function of crosswind 
distance (Problem 7). 

The values necessary to determine the isopleth 
half widths, y, are given in Table 7-3. 

Table 7-3 DETERMINATION OF ISOPLETH WIDTHS 
(PROBLEM 8) 



km 


ITy, 

m 


X (centerlinc), 
8 m~" 


X (isopleth) 


/,,. 


y, 
m 


X (centerline) 


0.5 


83 


3.8 X 10--"' 


0,263 


1.64 


136 


0.8 


129 


2.3 x 10- 1 


4,35 X 10~ a 


2.50 


323 


1.0 


157 


2.8x10-' 


3.53 X 10- 2 


2.59 


407 


2.0 


295 


1.4 x 10-' 1 


7.14 xlO- 2 


2.30 


679 


3.0 


425 


7.1 x 10- fl 


1.42X11T 1 


1.98 


842 


4.0 


540 


4.0 x 10-' 


0.250 


1.67 


902 


5.0 


670 


2.4 x 10" n 


0.417 


1.32 


884 


6.0 


780 


1.8 x 10~ r ' 


0.556 


1.08 


842 


7,0 


890 


i.4xio-< 


0.714 


0,82 


730 


8.0 


980 


1.1 x 10-' 


0.909 


0.44 


432 



The orientation of the x~axis will be toward 
225 close to the 'source, curving more toward 
210 to 215 azimuth at greater distances be- 
cause of the change of wind direction with 
height. The isopleth is shown in Figure 7-3. 

Since the isopleth approximates an ellipse, the 
area may be estimated by * ab where a is the 
semimajor axis and b is the semiminor axis. 



a 



86QO-350_ _ 4125 m 



2 



A (m 2 ) = TT (4125) (902) 

-=11.7xlOm 11 
or A = 11.7 km- 




SOUftCE 



Figure 7-3. Location of the 10" a g m" 3 ground-level iso- 
pleth (Problem 8). 

PROBLEM 9: For the conditions given in problem 
4, determine the profile of concentration with 
height from ground level to z = 450 meters at 
x = 1 km, y meters, and draw a graph of 
concentration against height above ground. 

SOLUTION: Eq, (3.1) is used to solve this prob- 
lem. The exponential involving y is equal to 1. 
At x = 1 km, tr y = 157 m, ^ = 110 m. (From 
problem 4) . 



Q 



151 



(110) 4 



3.5 x 10" g n 



Values for the estimation of x(z) are given in 
Table 7-4. 

PROBLEM 10: For the conditions given in prob- 
lem 4, determine the distance at which the 
ground-level centerline concentration equals the 
centerline concentration at 150 meters above 
ground. Verify by computation of x (x>0 3 0) 
and* (x,0,150). 

SOLUTION: The distance at which concentra- 
tions at the ground and at plume height are 
equal should occur where <r, = 0.91 H (bee 
Chapter 5). For B stability and H = 150 m, 
ff , = 0.91 (150) = 136 m occurs at x 1.2 km. 
At this distance a y = 181 m. 
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Table 7-4 DETERMINATION OF CONCENTRATIONS FOR 
VARIOUS HEIGHTS (PROBLEM 9) 



d. 



f. 



]z+H f 1 / z+H Vl 
"M ) j 



< 



01.36 


0.397 


1.36 


0.397 




0,794 


2 


.78 x 


ir j 


301.09 





.552 


1.64 


0.261 




0.813 


2.85 x 


10- 


60- 


-0.82 





.714 


1.91 


0.161 




0.875 


3.06 x 


10- 


90- 


-0.55 





.860 


2.13 


0.0929 




0.953 


3, 


.34 x 


10- 


1200,27 





.964 


2.45 


0.0497 


1.014 


3.55 x 


10- 


150 


0.0 


1 


.0 


2.73 


0,0241 




1.024 


3.58 x 


10r* 


180 


0.27 





.964 


3,00 


1.11 x 


10- 


0.975 


3,41 x 


1&~ 4 


210 


0.55 





.860 


3.27 


4.77 x 


10- 


0.865 


3.03 x 


10- 


240 


0.82 


0,714 


3.54 


1.90 x 


10-' 


0.716 


2.51 x 


10- 


270 


1.09 


0.552 


3.82 


6.78 x 


10- 


0.553 


1.94 x 


10-* 


300 


1.36 


0. 


397 


4.09 


2.33 x 


10- 


0.397 


1, 


39 x 


10-' 


330 


1,64 


0. 


261 


4.36 


7.45 x 


10- 


0.261 


9. 


14 x 


10- 


360 


1.91 


0. 


161 


4,64 


2.11 x 


10- 


0.161 


5.64 x 


10- 


390 


2.18 


0. 


0929 


4.91 


5.82 x 


10- 


0.093 


3. 


26 x 


io- 


420 


2.45 


0. 


0497 


5.18 


1.49 x 


10- 


0,050 


1. 


75 x 


10- 


450 


2.73 


0.0241 


5,45 


3.55x10- 


0.024 


8. 


40 x 


10- 



These values are plotted in Figure 7-4. 
500 
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Figure 74. Concentration as a function of height (Prob- 
lem 9). 



Verifying: 

x (x,0,0) 



Q 



TT (TV ov U 



exp 



151 f J^ 

181 (136) 4 exp [~~ 2 



1SQ 



4.88 x 10-* exp [ 
4.88 x 10- (0.546) 

2.7 x 10- g m- 




161 



2.181(136)4 
- 



-f exp 



lob 



F 1 1 

[ g-(2.21)| 



,,, 
-2.44 xIO- 1.0 - p 

_ 2A4 x ^ (1 _ Q + ^ x 1Q . 2) 

... v J 

= 2.44 x 10 - (1.087) 

_ 9 7 v 1 fl~4 tr m~U 

^./xiu gm 



PROBLEM 11: For the power plant in problem 4, 
what will the maximum ground-level concentra- 
tion be beneath the plume centerline and at 
what distance will it occur on a clear night with 
wind speed 4 m sec" 1 ? 

SOLUTION: A clear night with wind speed 4 in 
sec indicates E stability conditions. From Fig- 
ure 3-9, the maximum concentration should 
occur at a distance of 13 km, and the maximum 
XU/Qis 1.7x10- 



XHUIX 



Q 



1.7 x 10- x 151 



Q u ' 4 

= 6.4 x 10- g m- of S0 2 



PROBLEM 12: For the situation in problem 11, 
what would the fumigation concentration be the 
next morning at this point (x = 13 km) when 
superadiabatic lapse rates extend to include 
most of the plume and it is assumed that wind 
speed and direction remain unchanged? 

SOLUTION: The concentration during fumiga- 
tion conditions is given by Eq. (5.2) with the 
exponential involving y equal to 1. in this prob- 
lem. 

Q 



(x,0,0;H) = _, 



U 



For the stable conditions, which were assumed 
to be class E, at x = 13 km, a y = 520 m., and 
(7* = 90 m. Using Eq. (5.3) to solve for hj: 
hi = H + 2 a. -= 150 + 2 (90) = 330 m. 
From the horizontal spreading suggested by Eq. 
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ff , r = a y (stable) + H/8 = 520 + 19 = 539 

_ 151 

XP ~ ^ J (539) 330 

Note that the fumigation concentrations under 
these conditions are about 1.3 times the maxi- 
mum ground-level concentrations that occurred 
during the night (problem 11). 

PROBLEM 13: An air sampling station is located 
at an azimuth of 203 from a cement plant at a 
distance of 1500 meters. The cement plant re- 
leases fine particulates (less than 15 microns 
diameter) at the rate of 750 pounds per hour 
from a 30-meter stack. What is the contribution 
from the cement plant to the total suspended 
particulate concentration at the sampling sta- 
tion when the wind is from 30 at 3 m sec" 1 on 
a clear day in the late fall at 1600? 

SOLUTION: For this season and time of day the 
C class stability should apply. Since the sam- 
pling station is off the plume axis, the x and y 
distances can be calculated: 

x = 1500 cos T 1489 
y = 1500 sin 7 = 183 
The source strength is: 

Q = 750 Ib hr 1 x 0.126 fj? ' = 94.5 g sec" 1 
n Ib hr 1 

At this distance, 1489 m, for stability C, a y 
150 m, cr, 87. The contribution to the concen- 
tration can be calculated from Eq. (3.2): 



*r f\vr\ I _ 



exp 1 


lV t LLf ' 
7T Oy (T z U 

1 /H\"j 


-AJJ I 

[0 5 


2 \ 
/ 183 


2 U J J 

94 ' 5 c-p 


exp 


* 150 (87) 3 ' tp 
05 ( 3 Vl 


r 0.5 i 


V IbU 
fl 9.M 


' 5 \81 ) \ 

94.5 



1.23x10" 
exp [0.5 (0.345) 2 ] 

= 7,68x10-* (0.475) (0.943) 



PROBLEM 14: A proposed source is to emit 72 g 
sec- 1 of S0 a from a stack 30 meters high with 
a diameter of 1.5 meters. The effluent gases are 
emitted at a temperature of 250F (394K) 
with an exit velocity of 13 m sec"- 1 . Plot on log- 
log paper a graph of maximum ground-level 



concentration as a function of wind speeu or 
stability classes B and D. Determine the crit- 
ical wind speed for these stabilities, i.e., the 
wind speed that results in the highest concen- 
trations. Assume that the design atmospheric 
pressure is 970 mb and the design ambient air 
temperature is 20C (293K). 



SOLUTION: Using 
height equation: 



Holland's effective stack 



d 



u 

13 (1.6) 
u 

394 - 293 
394 

19.5 



1.5 + 2.68 x 10~ 3 p 



T a -T 



I 



u 

19.5 
u 

19.5 
u 

19.5 (2.5) 
u 



1.6 + 2.68 x 10- 3 (970) 

(1.5) 

1.5 + 2.6 (- 



101 
394 



1.5 



[1.5 + 2.6 (0.256) 1.5] 
[1.5 + 1.0] 



_ 48.8 
u 

The effective stack heights for various wind 
speeds and stabilities are summarized in Table 

7-5. 

Table 7-5 EFFECTIVE STACK HEIGHTS (PROBLEM 14) 



u, 
m sec- 1 


Class D 




Class B 




AH, h 
m 


i + AH, 
m 


1.35 AH r h + 

m 


1.15 AH, 
m 


0.5 


97.6 


127.6 


112.2 


142.2 


1.0 


18.8 


78.8 


56.1 


86,1 


1.5 


32,6 


62.6 


37.5 


67.5 


2 


24.4 


54.4 


28.1 


58.1 


3 


16.3 


46.3 


18.7 


48.7 


5 


9.8 


39.8 


11.3 


41.3 


7 


7.0 


37.0 


8.0 


38.0 


10 


4.9 


34.9 






20 


2.4 


32.4 







By use of the appropriate height, H, the maxi- 
mum concentration for each wind speed and 
stability can be determined by obtaining the 
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maximum *u/Q as a function of H and stability 
from Figure 3-9 and multiplying by the appro- 
priate Q/u. The computations are summarized 
in Table 7-6, and plotted in Figure 7-5. 



10" 



T T 




0.5 



1 2345 
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20 



Figure 7-5. Maximum concentration as a function of 
wind speed (Problem 14). 

Tabla 7-6 MAXIMUM CONCENTRATION AS A FUNCTION OF 
WIND SPEED (PROBLEM 14) 



Stability u, 
Class m sec~ 


H, 
1 m 


*"" 


Q/u,. 


E m 3 


B 0.5 


142,2 


8,0 x lO" 11 


144 


1.15x10-' 


1.0 


86.1 


2.0 x 1Q- E 


72 


1.44 xlO- 3 


1.5 


67.5 


3.1 x 10- fi 


48 


1.49xlQ- 3 <- 


2 


58.1 


4.1 x 10-' 


36 


1.48 x lO- 1 


3 


48,7 


5.7 x 10- 


24 


1.37 xlO- 3 


5 


41.3 


7.8xlO~ n 


14.4 


1.12 x 10-" 


7 


38.0 


8.7x10- 


10.3 


8.96 xlO- 1 


D 0.5 


127.6 


4.4 x 1Q-" 


144 


6.34 xlO- 1 


1.0 


78.8 


1.42x10"' 


72 


1.02 x 10- 3 


1,5 


62.6 


2,47x10-' 


48 


1.19 xlO- 3 


2 


54,4 


3.5x10"' 


36 


1.26xlO~ a <- 


3 


46.3 


5.1x10- 


24 


1.22 x 10-" 


5 


39.8 


7.3 x 10-' 


14,4 


1.05 x 10~ 3 


7 


37.0 


8.2xlO~ fi 


10.3 


8.45 x 10-* 


10 


34.9 


9.4x10-' 


7.2 


6.77x10-* 


20 


32.4 


1.1 xlO- 1 


3.6 


3.96x10-* 



The wind speeds that give the highest maximum 
concentrations for each stability are, from Fig- 
ure 7-5; B 1.5, D 2.0. 

PROBLEM 15: A proposed pulp processing plant 
is expected to emit % ton per day of hydrogen 
sulfide from a single stack. The company prop- 
erty extends a minimum of 1500 meters from 
the proposed location. The nearest receptor 



is a small town of 500 inhabitants 1700 meters 
northeast of the plant. Plant managers have 
decided that it is desirable to maintain 
concentrations below 20 ppb (parts per billion 
by volume), or approximately 2.9 x 10~ B g m" 3 , 
for any period greater than 30 minutes. Wind 
direction frequencies indicate that winds blow 
from the proposed location toward this town 
between 10 and 15 per cent of the time. "What 
height stack should be erected? It is assumed 
that a design wind speed of 2 m sec" 1 will he 
sufficient, since the effective stack rise will be 
quite great with winds less than 2 m sec" 1 . 
Other than this stipulation, assume that the 
physical stack height and effective stack height 
are the same, to incorporate a slight safety 
factor. 



SOLUTION: The source strength is: 

n _ 1000 Ib day" 1 x 453.6 gib- 1 
^ 86,400 sec day- 1 

FromEq. (4.2): 

0.117 Q 0.117 (5.25) 



= 5.25 g see" 1 



(2.9x10-") 2 



= 1.06 x 10* m 2 



At a design distance of 1500 meters (the limit 
of company property), <r y tr x = 1,06 x 10* gives 
a point from Figure 4-1 about 0.2 from Class C 
to Class D along the line x = 1500 m. From 
Figure 3-3, <r a = 80 for this stability. 
H = /2~ r s! = H3 meters 



PROBLEM 16: In problem 15 assume that the 
stack diameter is to be 8 ft, the temperature oi 
the effluent 250 F, and the stack gas velocity 
45 ft sec" 1 . From Holland's equation for effec- 
tive stack height and the method used in prob- 
lem 15, determine the physical stack height 
required to satisfy the conditions in problem 16. 
In estimating AH, use T n 68 F and p = 920 
mb. 

SOLUTION: First determine the relation between 
AH and u from Holland's equation. 

V H = 45 ft sec" 1 = 13.7 m sec" 1 

d = 8 ft 2.44 m 
T H 250F = 121C = 394K 
T ft = 68F - 20C = 293K 

p 920 mb 
v d 



u 

13.7 (2.44) 
u 

394-293 



1.5 + 2.68 x 10~ 3 p 



T 
AH 



1.5 + 2.68 x 10-" (920) 



394 



(2.44) 
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_ 33 - 4 [1.5 + (2.46) 0.256 (2.44)] 



60 sec min" 



u 

33.4 
u 

102 
u 



The relation between cr y o- z and u is: 

0.117 Q _ 0.117 (5.26) _ 
' ffz ~ y,i u 2.9 x 10~ 5 u 



2.12 x 1Q 4 
u 



The required computations using Figure 4-1 are 
summarized in Table 7-7: 

Table 7-7 REQUIRED PHYSICAL STACK HEIGHT AS A 
FUNCTION OF WIND SPEED (PROBLEM 16) 



u, AH, 
m sec" 1 m 


ffy 0. 


Stability to 
f Give o- y (r a at 
1500m 


m 


m 


L , 

H'-AH, 

m 


0.5 


204 


4.24 x 


10* 


0,9 


from 


A 


to 


B 


190 


269 


65 


1.0 


102 


2.12 x 


10* 


0.6 


from 


B 


to 


C 


120 


170 


68 


1.5 


68 


1.41 x 


10* 


0.9 


from 


B 


to 


C 


96 


136 


68 


2.0 


51 


1.06 x 


10* 


0.2 


from 


C 


to 


D 


76 


108 


57 


2.5 


41 


8.48 x 


10 n 


0.4 


from 


C 


to 


D 


64 


91 


50 


3.0 


34 


7.06 x 


10 a 


0.6 


from 


C 


to 


D 


56 


79 


45 


5.0 


20 


4.24 x 


10 3 




D 








42 


60 


40 


7.0 


15 


3.03 x 


10 3 


0.5 


from 


D 


to 


E 


34 


48 


33 


30.0 


10 


2.12 x 


10 3 




E 








28 


40 


30 


15.0 


7 


1.41 x 


10" 


0.5 


from 


E 


to 


F 


23 


33 


26 



The required physical height is 68 meters. 

PROBLEM 17: A dispersion study is being made 
over relatively open terrain with fluorescent 
particles whose size yields 1.8 x 10 10 particles 
per gram of tracer. Sampling is by membrane 
filters through which 9 x 10~ 8 m 3 of air is drawn 
each minute. A study involving a 1-hour release, 
which can be considered from ground-level, is to 
take place during conditions forecast to be 
slightly unstable with winds 5 m seer 1 . It is 
desirable to obtain a particle count of at least 
20 particles upon membrane filters located _ at 
ground-level 2.0 km from the plume centerline 
on the sampling arc 8 km from the source. What 
should the total release be, in grams, for this 
run? 

SOLUTION: The total dosage at the sampler is 
determined by the total sample in grama divided 
by the sampling rate: 



D T (gsecnr 3 ) = 



20 particles 



9 x 10~ 3 
1200 



m 3 



16.2 x 1Q 7 

D T 7.41 x 10- g sec nr a 
The total dosage is given in g sec m~ 3 from 

QT f 1 

; 0) ^ exp I FT 

17 U OV 0% I & 



D 



where Q- r is the total release in grams. 

mi f\ 7T U <Ty <T Z DT 

Therefore Q T = - 



\ J-f y Vl 

h 2 1; j 



For slightly unstable conditions (Class C) at 
x = 8 km, <r y = 690 m, a* = 310 m; y 2000 m, 

u = 5 m sec" 1 

TT 5 (690) 310 (7.41x10-) 

2000 
x 690 

24.9 



exp 




ri 



exp [0.5 (2.90) 2 ] 
24.9 



1.49 x 10~ a 
Q T = 1670 g 

No correction has been made for the facts that 
the release Is for 1 hour and the standard Devia- 
tions represent time periods of 3 to 15 minutes, 

PROBLEM 18: A release of 2 kg of fluorescent 
particles is made based on the results of the 
computation in problem 17. The conditions are 
class C stability and wind speed 5 rn sec *. The 
crosswind-integrated ground-level dosage along 
the 8-km arc is determined from the samplers 
along this arc to be 8.2 x 10' 1 g sec m" 2 . What 
is the effective * for this run? 



SOLUTION: 
given by: 



The crosswind-integrated dosage is 
2Q T 



a, u 



exp 






Since the source is at ground-level, the expo- 
nential has a value of 1. Solving for a t \ 



Y2TT DOWI 
2 (20QQ) 



1.8 x 10 10 particles g~ L 



4000 
10.28 

389m 



Example Problems 



PROBLEM 19: At a point directly downwind 
from a ground-level source the 3- to 15-minute 
concentration Is estimated to be 3.4 x 10~ 3 g 
m"'. "What would you estimate the 2-hour con- 
centration to be at this point, assuming no 
change in stability or wind velocity? 

SOLUTION: Using Eq. (5.12) and letting k = 3 
min, s = 2 hours, and p = 0.2; 



X 2 llullr 



120 y 

1 (0 


,4 x 10" 3 ) 

1 R v 


40 (d ' 
3.4 x 10-" 



2.09 



Letting k 15 min, s 2 hours, and p = 0,17 

X 8 Uonr = (-^-) "' T 3.4X10-" 

= 4^ (3.4xlO- 3 ) 



3.4xlQ- 3 
1.42 



2.4 x 



The 2-hour concentration is estimated to be 
between 1.6 x 10~' n and 2.4 x 10~ ;t g m~ ;t . 

PROBLEM 20: Two sources of S0 2 are shown as 
points A and B in Figure 7-6. On a sunny 
summer afternoon the surface wind is from 60 
at 6 m sec" 1 . Source A is a power plant emitting 
1450 g sec" 1 S0 a from two stacks whose physical 
height is 120 meters and whose AH, from Hol- 
land's equation, is AH (m) = 538 (m 2 sec^J/u 
(m sec" 1 )- Source B is a refinery emitting 126 g 
sec" 1 S0 2 from an effective height of 60 meters. 
The wind measured at 160 meters on a nearby 
TV tower is from 70 at 8.5 m sec" 1 . Assuming 
that the mean direction of travel of both plumes 
is 245, and there are no other sources of S0 2) 
what is the concentration of SO 2 at the receptor 
shown in the figure? 

SOLUTION: Calculate the effective height of 
Source A using the observed wind speed at 160 
meters. 



AH 



538 



63.3 



H A 120 + 63 -= 183 m 
QA = 1450 g sec~' 
H B = 60 m 
Q Il = 126gsec" 1 

For a sunny summer afternoon with wind speed 
6 m sec a , the stability class to be expected is C. 
The equation to be used is Eq. (3.2) ; 



SOURCE A 
x=!4.6 km 
y= 8,* 



RECEPTOR 




SCALE, km 
2 4 

Figure 7-6. Locations of sources and receptor (Problem 

20). 



x (x,y,0;H) 



Q 



exp 



TT try o- a U 

1 / H 



exp 



2 



For Source A, x = 24.6 km, y = 8.4 km 
ffy =_ 1810 m, <7 K = 1120 m, u = 8.5 m sec" 1 
1450 



XA 



exp I 0.5 




5.42 
exp [0.5 (0.164) 2 ] 

= 2.67x10") (2.11 xlO" 5 ) (0,987) 

For Source B, x = 13.0 km, y 4.0 km. 
o> 1050 m, CTK = 640 m, u = 7.0 m sec" 1 

_ 126 f n K f 400( > "\ ~ 

x " TT 1050 (640) 7 

f n ^ / 60 
exp I 0. 



- lJ 2 Ao' exp [-0,5 (3,81)^] 

exp [0.5 (0.0938) 2 ] 

= 8.5 x 10" (7.04 x 10"*) (0.996) 
XK 6,0 x 10" g m" a 
X XA + XB 0.56 x 10- + 6.0 x 10- 

6.6 x 10" g m ~ 3 
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PROBLEM 21: A stack 15 meters high emits 3 g 
sec- 1 of a particular air pollutant. The sur- 
rounding terrain is relatively flat except for a 
rounded hill about 3 km to the northeast whose 
crest extends 15 meters above the stack top. 
What is the highest 3- to 15-minute concentra- 
tion of this pollutant that can be expected on 
the facing slope of the hill on a clear night when 
the wind is blowing directly from the stack 
toward the hill at 4 m sec" 1 ? Assume that AH 
is less than 15 m. How much does the wind 
have to shift so that concentrations at this point 
drop below 10" 7 g rrr s ? 

SOLUTION: A clear night with 4 m sec' 1 indi- 
cates class E stability. Eq. (3.4) for ground- 
level concentrations from a ground-level source 
is most applicable (See Chapter 5). At 3 km 
for class E, <r y =- 140 m, <r a =- 43 m. 

Q 



" ~ TT 140 (43) 4 
X = 3.97 x 10- g nr 3 

To determine the crosswind distance from the 
plume centerline to produce a concentration of 
10~ 7 g m" 3 Eq. (3.8) is used: 



y- 2 In 



2 In 



(x,0,0) 



x (x,y,o) 

3.97 x 10" 



io- 7 

(21n397)^ 2 140 

(2x5.98)^140 

3.4G x 140 
= 484 m. 

484 



1/2 



(140) 



tan0 = 



3000 
= 9.2 



= 0.1614 



A wind shift of 9.2 is required to reduce the 
concentration to 10" 7 g m~ B . 

PROBLEM 22: An inventory of S0 2 emissions 
has been conducted in an urban area by square 
areas, 5000 ft (1524 meters) on a side. The 
emissions from one such area are estimated to 
be 6 g sec" 1 for the entire area. This square is 
composed of residences and a few small com- 
mercial establishments. What is the concentra- 
tion resulting from this area at the center of the 
adjacent square to the north when the wind as 
blowing from the south on a thinly overcast 
night with the wind at 2.5 m sec' 1 ? The average 
effective stack height of these sources is assumed 
to be 20 meters. 

SOLUTION: A thinly overcast night with wind 
speed 2.5 m sec- 1 indicates stability of class E. 



Example Problems 



(It may actually be more unstable, since this is 
in a built-up area. ) To allow for the area source, 
let (T,o = 1524/4.3 =- 354. For class E the vir- 
tual distance, x y 8.5 km. For x =- 1524 m, 
as = 28.5. For x + x y =10,024 m, a y =- 410 m. 

Q [ 1 /H 
^ - exp s- 



7T (Tj. 



6 




(28.5) 2.5~ 6XP 

= 6.:54xlO" 5 (0.783) 
x = 5.1x10" gmr :i 

PROBLEM 23: An estimate is required of the 
total hydrocarbon concentration 300 meters 
downwind of an expressway at 1730 on an over- 
cast day with wind speed 4 m see- 1 . The ex- 
pressway runs north-south and the wind is from 
the west. The measured traffic flow is 8000 
vehicles per hour during this rush hour, and the 
average speed of the vehicles is 40 miles per 
hour. At this speed the average vehicle is ex- 
pected to emit 2 x 10-' g sec- 1 of total hydro- 
carbons. 

SOLUTION: The expressway may be considered 
as a continuous infinite line source.^ To obtain 
a source strength q in grams sec 1 m l , the num- 
ber of vehicles per meter of highway must be 
calculated and multiplied by the emission per 
vehicle. 

Vehicles/meter = 

Flow (vehicles hour"') 

Average speed (miles hour 1 ) 1600 (mmile- 1 ) 



8QOQ 



= 1.25 x 1Q- 1 (vehicles irr l ) 



~~ 40x1600 

q 1.25 x 10-' (vehicles m" 1 ) x 2 x 1Q- 2 
(g sec" 1 vehicle" 1 ) 

Under overcast conditions with wind speed 4 m 
sec- 1 stability class D applies. Under D, at x 
300 meters, cr a = 12 m. From Eq. (5.18) : 

2q 



_ 2 (2.5x10"") 
~~ 2.507 (12) 4 
4.2 x 10" n g m~* of total hydrocarbons. 

PROBLEM 24: A line of burning agricultural 
waste can be considered a finite line source 150 
m long It is estimated that the total emission 
of orgamcs is at a rate of 90 g sec- 1 . What is the 
3- to 16-minute concentration of orgamcs at a 
distance of 400 m directly downwind from the 
center of the line when the wind is blowing at 
3 m sec-' perpendicular to the line? Assume 
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that it is 1600 on a sunny fall afternoon. What 
is the concentration directly downwind from one 
end of the source? 

SOLUTION: Late afternoon at this time of year 
implies slight insolation, which with 3 m sec"** 1 
winds yields stability class C. For C stability 
at x = 400 m, v y = 45 m, <r a = 26 m. 



Q 



90 



150 150 

Eq. (5.20) is appropriate. 

i n.r\\ ^q 



0.6 g sec 



X U,UU;U> : 
exp ( 0.5 p 2 ) dp 

D = y ~" 75 

P ' a y 45 

= 4-1.67 
X (400,0,0;0) 




45 




exp (0.5 p 2 ) dp 

= 6.14 xlO~ 3 (0.91) 
= 5.6xlO~ 3 gm- 3 

For a point downwind of one of the ends of the 
line: 



P, 0, p a = 



15Q 

45 




X (400,0,0:0) = 6.14 x 10- JI 

exp ( 0.5 p a ) dp 

= 6.14xlO- fl (0.4995) 
-=3.1xlO- 8 gnr It 

PROBLEM 25: A core melt-down of a power re- 
actor that has heen operating for over a year 
occurs at 0200, releasing 1.6 x 10" curies of 
activity (1 second after the accident) into the 
atmosphere of the containment vessel. This 
total activity can be expected to decay according 

to - -o- 2 . It is estimated that about 5.3 x 10* 



curies of this activity is due to iodine- 131, which 
has a half-life of 8.04 days, The reactor building 
is hemispherically shaped with a radius of 20 
meters. Assume the leak rate of the building is 
0.1% day' 1 . 



The accident has occurred on a relatively clear 
night with wind speed 2.5 m sec" 1 . What is the 
concentration in the air 3 kilometers directly 
downwind from the source at 0400 due to all 
radioactive material? due to iodine-131? 

SOLUTION: Source strength = leak rate x ac- 
tivity (corrected for decay) 

T i 4- 0.001 day- 1 

Leak rate Q/ .. nn 4 =r~ 

86400 sec day * 

-= 1.157 x 10- 8 sec' 1 
Source strength of all products 
Q A (curies sec' 1 ) = 1.157 x 10~ a (1.5 x 10) 
t (sec) 1 ~ ' 2 



be, (sec) J 
1.74 x 1Q- 2 



0.2 



To determine decay of materials with the half- 

life given, multiply by exp [ f 93t where t 
is time and L is half-life. \ ^ 

Source strength of I 181 . 



Qi (curies sec" 1 ) 
0.693 t 



1.157 x 10~ 8 (5.3 x 10*) exp 



L 
For I 181 L 6.95 x 10 fi sec 

0.693 t 



= 6.13 x 10-* exp 



6.95 x 10 n 



For a clear night with wind speed 2.5 m sec" 1 , 
class F applies. Approximate the spreading at 
the reactor shell by 2.15 o- y0 2.15 <r z0 = the 
radius of the shell -= 20 m o- y0 = <r z0 9-3 m. 
The virtual distances to account for this are: 
Xy = 250 m, X K = 560 m. 

At x 3000 m. x + x y 3250 m, <r r -= 100 m. 
x + X B =- 3560 m, <r a 29 m. 

Q Q 



X <x,0,0;0) 



TT (Ty 

4.4 x 10- G Q 



:U TT 100 (29) 2.5 



For concentration at 0400, 3000 m downwind 
due to all radioactivity, t = 7200 seconds. 

XA 4.4 x 10- (1.74 x 10~ a ) (7200)~- 2 
=- 7.66 x 10~ 7 (0.17) 

XA = 1.3 x 10~ 7 curies m" 3 

The concentration at 0400, 3000 m downwind 
due to I 131 is: 

Xi 4,4 x 10-" (6.13 x 10*) exp [0.997 x ID' 6 
(7200)] 
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=. 2.7 x 10~ K (1.0) The decay of I 131 is insig- 
nificant for 2 hours 

Xi = 2.7 x 10~ B curies m~ ;i 

PROBLEM 26: A spill estimated at 2.9 x 10" 
grams of unsymmetrical dimethyl hydrazine 
occurs at 0300 on a clear night while a rocket 
is being fueled. A circular area 60 meters in 
diameter built around the launch pad is revetted 
into squares 20 feet on a side to confine to as 
small an area as possible any spilled toxic liquids. 
In this spill only one such 20- by 20-foot area is 
involved. At the current wind speed of 2 m 
sec" 1 , it is estimated that the evaporation rate 
will be 1100 g sec" 1 . The wind direction is pre- 
dicted to be from 310 15 for the next hour. 
Table 7-8 gives the emergency tolerance limits 
for UDMH vapor. 



Table 7-8 EMERGENCY TOLERANCE LIMITS FOR 
VAPOR VERSUS EXPOSURE TIME 



Time, 
minutes 



Emergency Tolerance 
Limits, g m~ 3 



5 

15 
30 
60 



1.2 x 10-' 
8.6 x 10^- 
4.9 x 10- a 
2.5 x ID-" 



What area should be evacuated? 

SOLUTION: From Table 3-1, the stability class 
is determined to be Class F. This is not a point 
source but a small area source. Allowing 4.3 a y0 
to equal the width of the wetted area, 6.1 meters 
(20 feet), tfyit = 1.4 meters. In attempting to 
determine the virtual distance, x r , it is found to 
be less than 100 meters, and will be approxi- 
mated by 40 meters. The release will take: 



2.9 x 10 g 



1.1 x 10 :i gsec- ] 



2.64 x 10 ;t sec = 44 min. 



Therefore the concentration for an exposure 
time of 1 hour (2.5 x 10~ a g m~ :i ) is of main 
concern. 

The equation for calculation of downwind con- 
centrations is Eq, (3.4) : 



X (x,0,0;0) 

Of X 4- Xy, 



Q 



where o> is a function 



Values of the parameters and of x al ' e given in 
Table 7-9. 



Table 7-9 DETERMINATION OF CONCENTRATION AS A 
FUNCTION OF DISTANCE (PfiOBLEM 26) 



x, 
km 


a v 

in 


x + x v , 
Sim 


CTJ-, 
m 


X' 
g nr 3 


0,1 


2.3 


0.14 


5.5 


13.9 


0.3 


5.6 


0.34 


12.5 


2,5 


0.6 


9.7 


0.64 


22 


8.2 x ID' 1 


1 


14 


1.04 


35 


3.6 xlQ- 1 


3 


27 


3.04 


93 


7.0 x 10- 2 


6 


37 


6.04 


175 


2.7 x 10~ 2 


10 


47 


10.04 


275 


1.4 x io- a 



These values of x are graphed as a function of x 
in Figure 7-7. The downwind concentration 
drops below the critical value of 2.5 x 10~- at a 
distance of 6.5 km. 



100 



10 
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> 

SC 
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x 
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1 
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10 



Figure 7-7. Concentration of UDMH as a function of down- 
wind distance (Problem 26). 

Calculated widths within a given isopleth are 
summarized in Table 7-10. 

The maximum width of the area encompassed 
by an isopleth is about 140 meters from the 
downwind position. Since the wind direction is 
expected to be from 310itl5 , the sector at an 
azimuth of 115 to 145 plus a 140-meter rectan- 
gle on either side should be evacuated. 
See Figure 7-8, 



Example Problems 
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Table 7-10 DETERMINATION OF WIDTHS WITHIN 
ISOPLETHS (PROBLEM 26) 



x, 
km 


km 


cr v , 

m 


X (centerline), 
g jrr :t 


^(isoplelii) 


y 


y, 

in 


^ (centerline} 


IT r 


0.1 


0.14 


5.5 


13,9 


1.8 xHT 1 


3.55 


20 


0,5 


0,54 


19 


1.1 


2,27 x 10-- 


2.75 


52 


1.0 


1.04 


35 


3.6 x 10-' 


6.94 X 10- a 


2.31 


80 


2.0 


2,04 


66 


1.3x10-' 


1.92 x 10-> 


1.82 


120 


3.0 


3.04 


93 


7.0 x 10-* 


3.57 x 10~ ! 


1.44 


134 


4,0 


4.04 


120 


4.8 x 10~ a 


5.20 xlO-' 


1.14 


137 


5.0 


5.04 


149 


3.5 x 10-* 


7.14x10-' 


0.82 


122 


6.0 


6.04 


175 


2.7 x 10-' 


9.26 x 10- 1 


0.39 


68 




SCALE, km 



Figure 7-8. Possible positions of the 2.5 x 10~ a g m' 
isopleth and the evacuation area (Problem 26). 
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APPENDICES 



Appendix 1: ABBREVIATIONS AND SYMBOLS 
Abbreviations 

cal calorie 

g gram 

K degrees Kelvin 

ni meter 

rnb millibar 

sec second 

Syml)ol 

a ratio of horizontal eddy velocity to vertical 

eddy velocity 

c,, specific heat at constant pressure 
C y Sutton horizontal dispersion parameter 
C E Sutton vertical dispersion parameter 
d inside stack diameter at stack top 
DT (x,y,0;H) Total dosage 
e 2.7183, the base of natural logarithms 



f (0,S,N) 



h 
h, 

H 



K 

L 



frequency of wind direction for a given 
stability and wind speed class 
physical stack height 
height of the base of an inversion 
effective height of emission 
effective height of emission for a particular 
wind speed 

von Karman's constant, approximately equal 
to 0.4 

eddy diffusivity 

two uses: 1. the height of an air layer that is 
relatively stable compared to the 
layer beneath it; a lid 
2. the half-life of a radioactive 

material 

Sutton's exponent 
an index for wind speed class 
three uses: 1. Bosanquet's horizontal disper- 
sion parameter 

2. atmospheric pressure 

3. a dummy variable in the equa- 
tion for a Gaussian distribution. 

q two uses: 1. Bosanquet's vertical dispersion 

parameter 
2. emission rate per length of a line 

source 

emission rate of a source 
total emission during an entire release 
net rate of sensible heating of an air column 
by solar radiation 

the length of the edge of a square area source 
an index for stability 
a short time period 



t,,, 



T 
T, 
u 

Us 

v' 

v 

v* 

w' 

X 



x* 

y 

z 
z, t 
SO 
Sz" 

AH 

n 



"XL 



time required for the mixing layer to develop 

from the top of the stack to the top of the 

plume 

a time period 

ambient air temperature 

stack gas temperature at stack top 

wind speed 

a mean wind speed for the wind speed class N. 
horizontal eddy velocity 
stack gas velocity at the stack top 
a velocity used by Calder 
vertical eddy velocity 

distance downwind in the direction of the 
mean wind 

design distance, a particular downwind dis- 
tance used for design purposes 
the distance at which er = 0.47L 
a virtual distance so that a* (x x ) equals the ini- 
tial standard deviation, v xo 
a virtual distance so that o> (x r ) equals the ini- 
tial standard deviation, o> 
a virtual distance so that a K (x a ) equals the ini- 
tial standard deviation, <r 20 
crosswind distance 
height above ground level 
roughness parameter 

the rate of change of potential temperature 

with height 

the rise of the plume centerline above the stack 

top 

two uses: 1. wind direction azimuth or sector 

2. potential temperature 
3.1416 

ambient air density 

the standard deviation of azimuth (wind direc- 
tion) as determined from a wind vane or bi- 
directional vane 

the standard deviation of wind elevation angle 
as determined from a bi-directional vane 
the standard deviation in the downwind direc- 
tion of a puff concentration distribution 
an initial downwind standard deviation 
the standard deviation in the crosswind direc- 
tion of the plume concentration distribution 
an initial crosswind standard deviation 
the standard deviation in the vertical of the 
plume concentration distribution 
an effective ^ equal to 0.8 L 
an initial vertical standard deviation 
the vertical standard deviation of the plume 
concentration at a particular downwind dis- 
tance for the stability, S. 
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the angle between the wind direction and a 
line source 
concentration 
wi crosswind-integrated concentration 

a ground-level concentration for design pur- 
poses 

inversion break-up fumigation concentration 
concentration measured over a sampling time, 
t k 

(S maximum ground-level centerline concentra- 
tion with respect to downwind distance 



X K concentration measured over a sampling time 
t a 

y 

~ relative concentration 
y 

relative concentration normalized for wind 
" speed 

X (x,y,z;H) concentration at the point (x, y, z) 
from an elevated source with effective 
height, H. 

X (x,o) the long-term average concentration at 
distance x, for a direction from a source. 
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Appendix 2: CHARACTERISTICS OF THE 
GAUSSIAN DISTRIBUTION 

The Gaussian or normal distribution can be de- 
picted by the bellshaped curve shown in Figure A-l. 
The equation for the ordinate value of this curve is: 



exp 



(A.1) 



Figure A-2 gives the ordinate value at any distance 
from the center of the distribution (which occurs 
at x) This information is also given in Table A-l. 
Figure A-3 gives the area under the Gaussian curve 
from ^ to a particular value of p where p = 



This area is found from Eq. (A.2): 



Area ( x to p) = 



exp ( 0.5 p' J ) dp 




(A.2) 



Figure A-4 gives the area under the Gaussian 
curve from p to -f p. This can be found from Eq. 
(A.3): 



Area ( p to -}-p) 
exp ( 0.5 p 2 ) dp 




(A.3) 




-3 



Figure A-l. The Gaussian distribution curve. 
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Figure A-2. Ordinate values of the Gaussian distribution. 
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Figure A-3. Area under the Gaussian distribution curve from to p. 



Appendix 




0.01 



0.5 1 2 5 10 20 30 40 50 60 70 80 90 95 98 99 

+p 

=- axp (-0.5 p 2 ) dp 

P 



99.8 99.99 



Figure A-4. Area under the Gaussian distribution curve between p and +p. 
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Appendix 3: SOLUTIONS TO EXPONENTIALS 

Expressions of the form exp [-0.5 A 8 ] where 
A is H/a, or y/tfv frequently must be evaluated. 
Table A-T'JvesB as a function of A where B - p 
lo.5 A"]- The sign and digits to the right of the 
E are to be considered as an exponent o 10. For 
example, if A is 3.51, B is given as 2.11E - Od 
which means 2.11 x 10 * 
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Appendix 4: CONSTANTS, CONVERSION 
EQUATIONS, CONVERSION TABLES 

Constants 

e = 2.7183 -J_ = 0.3679 



= 3.1416 L- = 0.3183 



TT 



27r = 6.2832 i = 0.1592 

2?r 

57= 2.5066 -= = 0.3989 

V27T 

^r = 0.7979 



(Sir) 3 /'-= 15.75 

Conversion Equations and Tables 

T(C) =5/9 <T(F) 32) 
T(K) =T(C) +273.16 
T(F) (9/6T(C) ) +32 
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